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ABSTRACT 



Present low energy neutrino oscillation data are elegantly explained by neutrino oscil- 
lation hypothesis with very small masses (< leV) of the light neutrinos. These masses 
can be either Dirac or Majorana. Small Majorana masses of the light neutrinos, how- 
ever, can be generated through the seesaw mechanism without any fine tuning. This can 
be achieved by introducing right handed neutrinos into the electroweak model which are 
invariant under all gauge transformations. The Majorana masses of these right handed 
neutrinos are free parameters of the model and are expected to be either at TeV scale 
or at a higher scale. This indicates the existence of new physics beyond Standard Model 
(SM) at some predictable high energy scale. 

Beyond SM baryogenesis via leptogenesis is an attractive scenario that links the 
physics of right handed neutrino sector with the low energy neutrino data. Majorana 
mass of the neutrino violates lepton (L) number and thus provides a natural path to 
generate L-asymmetry. The leptogenesis occurs via the out of equilibrium decay of heavy 
right handed Majorana neutrinos to SM leptons and Higgses. Assuming a normal mass 
hierarchy in the right handed heavy Majorana neutrino sector the final L-asymmetry is 
given by the CP- violating decays of lightest right handed Majorana neutrino, N±. A 
partial L-asymmetry is then transformed to the baryon (B) asymmetry via the non per- 
turbative sphaleron processes. 

We divide the thesis into two parts. In part-I, we study baryogenesis via leptogenesis 
in a thermal scenario, while part-II of the thesis is devoted to study the same in a non- 
thermal scenario. In both scenarios we discuss bounds on the mass scale of right handed 
heavy Majorana neutrinos from the leptogenesis constraint. Moreover, we divide the 
phenomenological models into two categories, type-I and type-II, depending on the seesaw 
mechanism used to generate the light Majorana neutrino masses. 

Part-I of the thesis begins with a brief introduction to type-I seesaw models. In this 
model, the Majorana mass matrix of the light neutrinos is given by m u = moM^m^, 
where ma is the Dirac mass matrix of the neutrinos. On the other hand, in the type- 
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II seesaw models the Majorana mass matrix of the light neutrinos is given by m v = 
Ml — mjyM^m^), where the additional mass M L , in contrast to type-I case, is provided 
by the vacuum expectation value of the triplet A^. The two terms, Ml and m^M^m^, 
contributing the neutrino mass matrix m v are called type-II and type-I respectively. Irre- 
spective of the magnitudes of type-II and Type-I terms, it is shown that in a hierarchical 
mass basis of right handed Majorana neutrinos the leptogenesis constraint gives rise a 
lower bound on the mass scale of N\ to be M\ > O(10 8 ) GeV, assuming that the CP- 
violating decays of Aq produces the observed P-asymmetry via the leptogenesis route. 
Numerically we check the compatibility of this bound with the low energy neutrino oscil- 
lation data. 

As a specific example of type-II seesaw models, we consider Left-Right symmetric 
model with spontaneous CP- violation. The Lagrangian of this model is CP-invariant 
and the Yukawa couplings are real. Due to spontaneous breaking of the gauge symmetry, 
some of the neutral Higgses acquire complex vacuum expectation values, which lead to 
CP-violation. In this model, we identify the neutrino Dirac mass matrix with that of 
charged lepton mass matrix. We assume a hierarchical spectrum of the right handed 
neutrino masses and derive a bound on this hierarchy by assuming that the decays of Aq 
produces the observed P-asymmetry via the leptogenesis route. It is shown that the mass 
hierarchy we obtain is compatible with the current neutrino oscillation data. 

The bound on the mass scale of Aq, in production of L-asymmetry through it's CP- 
violating decays to SM Higgs and lepton, is > 10 8 GeV which is far above the current 
accelerator energy range and beyond the reach of the next generation accelerators. How- 
ever, these scenarios are well motivated by the current status of low energy neutrino 
oscillation data. An elegant alternative is to consider mechanisms which work at TeV 
scale and be consistent with the low energy neutrino data. We assume that the required 
asymmetry of the present Universe is produced during the B—L gauge symmetry breaking 
phase transition. Below Tb_l, the scale of B — L symmetry breaking phase transition, the 
preservation of L-asymmetry constrains the mass scale of Aq, to be O(10)TeV if the Dirac 
masses of the light neutrinos are of 10 2 order smaller than the charged lepton masses. By 
solving the required Boltzmann equations we check the compatibility of lOTeV scale right 
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handed neutrino with the low energy neutrino oscillation data. We discuss a scenario for 
the production of large L-asymmetry during the B — L gauge symmetry breaking phase 
transition. 

In part-II of the thesis, we discuss soliton-fermion systems in gauge theories. Solitons 
emerge as the time independent solutions of non linear wave equations in classical gauge 
theories. However, their interactions with fermions lead to a curious phenomenon of 
fractional fermion number. We have considered the possibility of fermion fractionization 
in various toy models and its implication for stabilizing otherwise metastable solitons. 

A typical solitonic solution in 3+1 dimensional gauge theory is cosmic string. It is 
a 1 + 1 dimensional extended object. During the early Universe phase transitions such 
objects are formed as topological defects. These objects are highly non-thermal and carry 
a fraction of energy of the Universe in their core called false vacuum. The decay of these 
objects produces quanta of massive particles, which may survive for long times and hence 
can provide a link between the early Universe and recent cosmology. In particular, we 
study baryogenesis via the route of leptogenesis. 

We study the contribution to the baryon asymmetry of the Universe (BAU) due to 
decay of heavy right handed Majorana neutrinos released from closed loops of B — L 
cosmic strings in the light of current ideas on light neutrino masses and mixings implied 
by atmospheric and solar neutrino measurements. We have estimated the contribution to 
BAU from cosmic string loops which disappear through the process of (a) slow shrinkage 
due to energy loss through gravitational radiation — which we call slow death (SD), 
and (b) repeated self-intersections — which we call quick death (QD). We find that for 
reasonable values of the relevant parameters, the SD process dominates over the QD 
process as far as their contribution to BAU is concerned. It is shown that for the B — L 
symmetry breaking scale ^b-l ~ 1-7 x 10 11 GeV the SD process of cosmic string loops 
contribute significantly to the present BAU. 

Keywords: seesaw mechanism,type-I seesaw models, type-II seesaw models, CP- 
asymmetry, leptogenesis, baryogenesis, soliton, cosmic strings, fermion zero modes. 
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Chapter 1 
Introduction 



One of the challenging problems in theoretical physics is baryogenesis. Since we live in a 
baryonic Universe it is worth investigating. Moreover, baryogenesis plays an important 
role in the interface of particle physics and cosmology and thus provides a scope to link 
them. Though there are a few evidences regarding the presence of antibaryons, they are 
tenuous. In particular, the presence of antiproton in the cosmic ray shower is one in 
O(10 4 ). 

The present Universe is electrically neutral. This is an indication of the U(l) em sym- 
metry of the present Universe. Therefore, without loss of generality we assume it holds 
since the Big-Bang. If baryon (B) number and lepton (L) number were absolutely con- 
served by all possible interactions occurring in the early Universe, then the total B and 
L numbers of the present Universe must simply reflect their apparently arbitrarily im- 
posed initial values. A plausible guess would be that the initial B and L numbers were 
exactly zero. That means each particle has its own antiparticle carrying an equal and op- 
posite quantum number thus maintaining a charge neutrality since the Big-Bang. Then 
the questions arise "where are the antibaryons (antimatter) ?" and " why the present 
Universe is baryonic (matter) ?". 

Assuming a highly symmetric state in the early Universe, a matter-antimatter asym- 
metry can be dynamically generated in an expanding Universe if the particle interactions 
and the cosmological evolution satisfy Sakharov criteria |T], i.e. 

• Baryon number violation. 
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• C (charge conjugation) and CP (charge conjugation plus parity)-violation. 

• Out of thermal equilibrium. 

Based on these criteria, several mechanisms have been put forwarded since late seven- 
ties. One of the early proposals is Grand Unified Theory (GLT)-baryogenesis [21 IS]- Since 
B — L is a gauge symmetry of most of the known GUT models, any asymmetry produced 
at the GUT scale will be erased. Thus the GUT solution of baryogenesis is unlikely to 
be true. During the early Universe phase transitions the last opportunity to produce the 
baryon asymmetry is the electroweak {EW) phase transition. The strategy is to assume a 
first order phase transition to ensure an epoch of non-equilibrium evolution jHEHHI, during 
which the B, C and CP violating effects must take place, satisfying the Sakharov crite- 
ria pp. However, the thermodynamics of EW phase transition indicates that a first order 
phase transition is unlikely [7] thus making baryogenesis unfeasible within the Standard 
Model (SM). 

In the thermal era of the early Universe a plausible explanation of the observed B- 
asymmetry of the Universe (BAU) is that it arose from a L-asymmetry [HJ E3 ITT] . 
The conversion of the L-asymmetry to the 5-asymmetry then occurs via the high tem- 
perature behavior of the B + L anomaly of the SM [12] IT3] ITl] . This is an appealing 
route for several reasons. First, the extremely small neutrino masses, suggested by the 
solar and atmospheric ^H] neutrino anomalies and the KamLAND experiment [T7] . 
point to the possibility of Majorana masses for the neutrinos generated by the seesaw 
mechanism [IE] that involves the right handed heavy Majorana neutrinos. This suggests 
the existence of new physics at a predictable high energy scale. Since the Majorana mass 
terms violate lepton number they can generate L-asymmetry in a natural way. Second, 
most particle physics models incorporating the above possibility demand new Yukawa 
couplings and also possibly scalar self-couplings; these are the kind of couplings which, 
unlike gauge couplings, can naturally accommodate adequate CP violation, one of the 
necessary ingredients pQ for generating the BAU. 

Most proposals along these lines rely on out-of-equilibrium decay of the heavy Ma- 
jorana neutrinos to generate the L-asymmetry. In the simplest scenario a right handed 
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neutrino per generation is added to the SM (HI E]. They are coupled to left handed 
neutrinos via Dirac charged lepton mass matrix JS]- Since the right handed neutrino is 
a singlet under SM gauge group a Majorana mass term (Mr) can be added to the La- 
grangian. Diagonalization of neutrino mass matrix leads to two Majorana neutrino states 
per generation: a light neutrino state (mass ~ m 2 D /M^) which is almost left handed and 
a heavy neutrino state (mass ~ Mr) which is almost right handed. This is called type-I 
seesaw mechanism JTHj in which there are no Majorana mass terms for the left handed 
fields. The wide class of models in which the light neutrinos derive their masses via type- 
I seesaw mechanism are called type-I seesaw models. In such models the right handed 
neutrinos do not possess any gauge interaction. An appealing way to solve this problem 
is to extend the SM by the inclusion of an extra ?7(1)b-l gauge symmetry [TT] . At a 
high scale, the singlet Higgs field x acquires a vacuum expectation value (VEV) and thus 
breaking the B — L gauge symmetry. The VEV of x provides heavy Majorana masses to 
the singlet right handed neutrinos through the Yukawa couplings. However, the B — L 
gauge symmetry in such models is quite ad hoc. 

An alternative is to consider models in which the B — L gauge symmetry emerges 
naturally. One of the possibilities is the Left-Right symmetric model SU(2) l x SU(2)r x 
U(1)b-l- ^ can be embedded in the higher gauge groups of most of the known Grand 
Unified Theories (GUTs) [19 . In such models Majorana masses, Ml, for left handed 
neutrinos occur in general, through the VEV of the triplet Al ffiH Effi 1^1 1^4] . The 
diagonalization of neutrino mass matrix in such models gives also a light and a heavy 
neutrino state per generation. The heavy neutrino state has mass ~ Mr but the light 
neutrino mass is ~ (Ml — m 2 D /MR). The two contributions to the light neutrino mass 
matrix, m 2 D /Mn and Ml are called type-I and type-II terms respectively. Such models 
in which both type-I and type-II terms contributing the light neutrino mass matrix are 
called type-II seesaw models. Because of B — L is a gauge charge of such models, no 
primordial B — L can exist. Further, the rapid violation of the B + L conservation by the 
anomaly due to the high temperature sphaleron fields erases any B + L generated earlier. 
Thus the L-asymmetry must have been produced entirely during or after the B — L gauge 
symmetry breaking phase transition. 
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The goal of the present neutrino oscillation experiments is to determine the nine pa- 
rameters in the leptonic mixing matrix assuming that the neutrino masses are Majorana 
in nature. The set of parameters include three light neutrino masses, three mixing angles 
and three phases which include one Dirac and two Majorana. At present the neutrino 
oscillation experiments able to measure the two mass square differences, the solar and 
the atmospheric, and three mixing angles with varying degrees of precision, while there 
is no information about the phases. The Majorana phases can be investigated in neutri- 
noless double beta decay experiments, while the Dirac phase can be investigated in the 
long base line neutrino oscillation experiments. Moreover, it is difficult to constrain the 
parameters of the right handed neutrinos from the low energy neutrino data. However, 
several attempts |2Sj have been made by inverting the neutrino mass matrix in type-I 
seesaw models. 

Baryogenesis via leptogenesis provides an attractive scenario to link the physics of 
right handed neutrino sector with the low energy neutrino data j2E]- We assume that 
the mass basis of right handed Majorana neutrinos, Ni, i — 1,2, 3, is diagonal. In this 
basis, we further assume that the mass spectrum of right handed Majorana neutrinos is 
in normal hierarchy, Mi < M 2 < M 3 . In this scenario, while the heavier neutrinos N 2 and 
N 3 decay, the lightest right handed Majorana neutrino Ni is still in thermal equilibrium. 
Thus any L-asymmetry produced by the decay of N 2 and N 3 is erased by the lepton 
number violating interactions mediated by N\. Hence the final L-asymmetry is given by 
the CP- violating decays of N\. The required L-asymmetry constrains the mass scale of 
Nt, in the type-I seesaw models, to be M x > 10 9 GeV [22 123 123 EH- On the other 
hand, in type-II seesaw models, where the Majorana mass term for left handed fields 
also contribute to the neutrino mass matrix, this bound can be reduced by an order of 
magnitude (22 E2]- This bound is well compatible with the low energy neutrino oscillation 
data. 

As a specific example of type-II seesaw models, we study Left-Right symmetric model. 
In this model, we (33] consider a special case in which CP-asymmetry arises through 
the spontaneous symmetry breaking [34] . The Lagrangian of the model is CP-invariant 
which demands that all the Yukawa couplings should be real. In this scenario the vacuum 
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expectation values (VEVs) of the neutral Higgses are complex which lead to CP- violation. 
In the Left-Right symmetric model, there are four complex neutral scalars which acquire 
VEVs. However, the unbroken global U(l) symmetries associated with SU(2)l and 
SU{2)r gauge groups allow two of the phases to be set to zero. Using the remnant U(l) 
symmetry after the breaking of S , f/(2)^, one phase choice is made to make the VEV of 
and hence the mass matrix of right handed neutrinos, real. The phase associated with 
the other U(l) symmetry can be chosen to achieve two different types of simplification of 
neutrino mass matrix. In the type- II choice, the ml is made real leaving the CP- violating 
phase purely with ml 1 . In this phase convention, we derive a lower bound on the mass 
scale of N\ from the leptogenesis constraint by assuming a normal mass hierarchy in 
the right handed neutrino sector. It is shown that the mass scale of N\ satisfy the 
constraint M l > 10 8 GeV [33], which is in good agreement with the lower bound on Mi 
in generic type-II seesaw models |31| I32j. In the type- 1 phase choice, only the type-I term 
contains CP-violating phase leaving type-II term real. This allows us to derive an upper 
bound on the heavy neutrino mass hierarchy from the leptogenesis constraint. In order to 
achieve the observed baryon asymmetry of the present Universe, it is found that the mass 
hierarchy of right handed neutrinos satisfy the constraint M 2 /Mi < 17 and M 3 /Mi < 289 
simultaneously [33] ■ Numerically we verified that these bounds are compatible with the 
low energy neutrino oscillation data for all values of Mi > 10 8 GeU as implied by the 
lower bound on Mi in type-II phase convention. 

The bound on Mi, in production of L-asymmetry through the CP- violating decays of 
thermally generated Ni, is > 10 8 GeV which is far above the current accelerator energy 
range and beyond the reach of the next generation accelerators. However, these scenarios 
are well motivated by the current status of low energy neutrino oscillation data. An 
alternative is to consider mechanisms which work at the TeV scale and may rely on the 
new particle content implied in supersymmetric extensions of the SM [33]. The Minimal 
Supersymmetric SM (MSSM) holds only marginal possibilities for baryogenesis. The 
Next to Minimal or NMSSM possesses robust mechanism for baryogenesis [3B] however 
the model has unresolved issues vis a vis the [i problem due to domain walls [37]. However 
its restricted version, the nMSSM is reported HH EH EH EJ 112] to tackle all of the 
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concerned issues. 

It is worth investigating other possibilities, whether or not supersymmetry is essential 
to the mechanism. We therefore studied jlHj the consequence of heavy Majorana neutrinos 
given the current knowledge of light neutrinos. The starting point is the observation [H| 
143] that the heavy neutrinos participate in the erasure of any pre-existing asymmetry 
through scattering as well as decay and inverse decay processes. Estimates using general 
behavior of the thermal rates lead to a conclusion that there is an upper bound on the 
temperature Tb_l a t which B — L asymmetry could have been created. This bound is 
Tb-l ~ !0 13 GeV x (leV/m u ) 2 , where m u is the typical light neutrino mass. We extend 
this analysis by numerical solution of the Boltzmann equations and obtain regions of 
viability in the parameter space spanned by m\-Mi, where fhi is called effective neutrino 
mass parameter. We find that our results are in consonance with [33] where it was argued 
that scattering processes provide a weaker constraint than the decay processes. If the 
scatterings become the main source of erasure of the primordial asymmetry then the 
constraint can be interpreted to imply Tb-l < Mi- Further, this temperature can be as 
low as TeV range with fh\ within the range expected from neutrino observations. This 
is compatible with seesaw mechanism if the "pivot" mass scale is two order smaller than 
that of the charged leptons. We conjecture that the hypothesis of TeV scale right handed 
neutrino can be verified in near future and hence an indirect evidence of generating the 
baryon asymmetry at TeV scale |43j . 

Part-II of this manuscript is devoted to study the formation and evolution of topolog- 
ical defects |4l>j, i n particular cosmic strings, in the early Universe. The corresponding 
consequences have been studied in details along the direction of baryogenesis via the route 
of leptogenesis. 

Topological defects arise as the solitonic solutions in gauge theories. 'Solitons' or 
'solitary waves' are the time independent solutions of non-linear wave equations in classical 
field theories. The prime among them is A0 4 theory. In 1+1-dimensions the solitonic 
solutions in A0 4 theory are called 'kinks'. Note that, kink is purely a solitary wave 
but not a soliton. However, in field theory the distinction between them is completely 
blurred. Each solitonic solution is designated by a number called the 'topological charge' 
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or 'winding number'. The topological charge is nothing but the boundary conditions 
imposed on the field which is conventionally different from the Noether charge comes 
from the continuous global symmetry associated with the theory. 

In Quantum Field Theory (QFT), solutions of Dirac equation in the presence of 
solitonic objects lead to a curious phenomenon of 'fractional fermion number' [37[|. This is 
because of the existence of degenerate zero energy modes of fermions while quantized in the 
background of a solitonic vacuum. In contrast to it, in the translational invariant vacuum 
there is no zero energy solutions of Dirac equation and therefore fermions are quantized by 
integral unit. The fractionally charged solitonic states are therefore superselected from 
the normal vacuum and are not allowed to decay in isolation. Note that the fermion 
number we are talking about is the eigenvalue of the number operator in QFT. 

An inevitable feature of the early Universe phase transitions is the formation of 
topological defects [IB). In particular we shall deal with cosmic strings. The break- 
down of any U(l) gauge symmetry to Z ensures the formation cosmic strings since 
TV i (Z7(l)/Z) = 7r (Z) = Z. These defects are extended objects and are not distributed 
thermally. Therefore, the decay of these objects can be a non-thermal source of massive 
particles. Moreover, the cosmic strings formed at a phase transition can also influence 
the nature of a subsequent phase transition that may have important implications for the 
generation of BAU 1381133] . 

An important feature of these cosmic strings is that during their formation they trap 
zero modes of fermions [HO] which are well predicted [HI] . These fermionic zero modes 
induce fractional fermion number (\n\/2), where n is the winding number of a string. If 
n is odd then the induced fermion number on the string is half-integral. Therefore, it 
is superselected from the translation invariant vacuum where the eigenvalues of the 
number operator possesses integral fermion number. So a string of half integer fermion 
number can not decay in isolation. However, this conclusion may not be true for closed 
loops which are chopped off from the infinite straight strings. This remains an important 
open question for cosmology. 

There exist both analytical as well as numerical studies of the evolution of cosmic 
string networks in the early Universe. These suggest that the string network quickly 
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enters a scaling regime in which the energy density of the strings scales as a fixed fraction 
of the energy density of radiation in the radiation dominated epoch or the energy density 
of matter in the matter dominated epoch. In both cases the energy density scales as the 
inverse square of cosmic time t. In this regime one of the fundamental physical process 
that maintains the strings network to be in that configuration is the formation of sub- 
horizon size closed loops which are pinched off from the network whenever a string segment 
curves over into a loop, intersecting itself. 

In many scenarios |53J EH EH1 EE1 EZj it has been studied that the decaying, collapsing, 
or repeatedly self-intersecting closed loops of such cosmic strings provide a non-thermal 
source of massive particles that "constitute" the string. The decay of these massive 
particles give rise to the observed B-asymmetry or at least can give a significant con- 
tribution to it. Assuming that the final demise of each string loop produces 0(1) right 
handed neutrino, N\, the observed 5-asymmetry [SH] requires the constraint on its mass 

1 /2 

Mi < 2.4 x 10 12 (?7b-l/10 13 GeV) GeV, t] B -l being the scale of £/(1)b-l gauge symme- 
try breaking phase transition. Here we have assumed that at a temperature above the 
mass scale of N\ there is no lepton asymmetry. A net asymmetry has been produced 
just below the mass scale of N± by it's CP-violating decay to SM Higgs and lepton. In 
order to take into account the wash out effects we solve the required Boltzmann equa- 
tions [5J El EH EH] by including the right handed neutrinos of cosmic string origin as 
well of thermal origin. It is shown that the delayed decay of cosmic string loops over 
produce the baryon asymmetry of the Universe for 7/b-l > O(10 n ) GeV jHUj. This, On 
the other hand, gives an upper bound on the CP-violating phase to produce the required 
asymmetry of the present Universe. 

The rest of the manuscript is organized as follows. In chapter 2, we briefly discuss the 
thermal baryogenesis via the route of leptogenesis in type-I seesaw models. In chapter 3, 
thermal baryogenesis via the route of leptogenesis in type-II seesaw models is discussed 
in detail. In both cases we find that the scale of operation of leptogenesis should be very 
high, of the order 0(> 10 8 ) GeV. In chapter 4, we propose the possibility that TeV scale 
masses for the right handed heavy neutrinos are consistent with the seesaw and may yet 
suffice to explain the baryon asymmetry. A relevant model is also summarized. The second 
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part of the thesis is devoted to study the formation and evolution of topological defects 
in the early Universe. In chapter 5, we briefly introduce the soliton-fermion systems in 
QFT. We than discuss the consequences of quantization of fermions in the background 
of solitons. The same hypothesis is extended to the case of cosmic strings in chapter 6. 
In chapter 7, evolution of B — L cosmic strings is discussed in greater detail. In a scaling 
regime the decay of the massive particles emitted from the cosmic string loops produce 
a baryon asymmetry via the leptogenesis route. The consequences for the energy scale 
of leptogenesis is discussed. The thesis ends by a summary and conclusions in chapter 8. 
Through out this manuscript we use natural units and set c — h — ks — 1- 
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Part-I 

Baryogenesis via Leptognesis 

in 

Thermal Scenario 
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Chapter 2 



Thermal leptogenesis in type-I 
seesaw models and bounds on 
neutrino masses 

2.1 Heavy Majorana neutrinos: The physics beyond 
standard model 

Within the SM the left handed neutrinos are massless. However, the present evidence 
from the neutrino oscillation experiments DZ3 El suggests that the left handed neutri- 
nos possess small masses (< leV). These masses can be either Dirac or Majorana. Small 
Majorana masses of the light neutrinos, however, can be generated via seesaw mecha- 
nism [IB] that involves the right handed heavy Majorana neutrinos. This indicates the 
existence of new physics beyond SM at a predictable high energy scale. 

In the simplest scenario a massive right handed neutrino iV^ of mass Mr per generation 
is added to the SM. They are coupled to the left handed neutrinos via Dirac mass matrix 
m D . Since the right handed neutrino is a singlet under SM gauge group, N R —N R coupling 
can be added to the SM Lagrangian. This gives a Majorana mass to the right handed 
neutrino. On the other hand, the vl — coupling is not allowed by the SM Lagrangian 
as it violates the lepton number by two units. Therefore, the upper left 3x3 block of the 
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6x6 neutrino mass matrix is zero; see for example section 2.2. The diagonalization of the 
neutrino mass matrix thus leads to two Majorana neutrino states per generation: a light 
neutrino state of mass ~ m D /Mn which is almost left handed and a heavy neutrino state 
of mass ~ Mr which is almost right handed. This is called type-I seesaw mechanism. 
The class of models in which the Majorana masses of the light neutrinos are obtained via 
this mechanism are called type-I seesaw models. 



2.2 Type-I seesaw mechanism and neutrino masses 

To generate the light neutrino masses via type-I seesaw mechanism we add a massive right 
handed neutrino per generation to the SM Lagrangian. For three generation of neutrinos, 
the terms in the Lagrangian for the massive right handed Majorana neutrinos are taken 
to be 



C Nr = iN Ri YduN Rl + 



v L im Dij N R3 + ^N c Ri M RlJ N Rj + H.C. 



(2.1) 



where i, j = 1, 2, 3 are family indices. Without loss of generality we can choose a basis in 
which the right handed Majorana neutrinos are diagonal. In this basis the neutrino mass 
matrix in 6 x 6 block is given as 

mo 



ml M f , 



(2.2) 

I'd 1v1 r' 

Diagonalizing the neutrino mass matrix (j2.2j) into 3x3 blocks we get the mass matrix 
for the light neutrinos to be 

1 rp 

m v = -m D ——m D 

= -v 2 h-^-h T , (2.3) 

Mr K J 

where v is the vacuum expectation value of SM Higgs and h is the relevant Yukawa cou- 
pling matrix. From equation ()2.3|) . it is clear that larger the mass scale of the right handed 
neutrinos the smaller is the mass scale of the left handed neutrinos. Diagonalization of 
the light neutrino mass matrix m u , through the lepton flavor mixing matrix Upmns |61j . 
gives us three light Majorana neutrino states. The eigenvalues are 

U]m v Ul = -dia(mi,m 2 ,m 3 ) = -D m , (2.4) 
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where m 1; m 2 , m 3 are chosen to be real. Combining equations (|2.3|) and (J2.4j) we get the 
diagonal mass matrix 

Dm = v 2 Ulh^-h T Ul. (2.5) 

The smallness of the light neutrino masses imply that the mass scales of right handed 
neutrinos, Ni,N 2 and N 3 exist at a high scale and beyond the energy range of current 
accelerators. However, their decay to SM particles may have significant consequences. 
In particular, we consider baryogenesis via the route of leptogenesis through the out-of- 
equilibrium decay of the right handed heavy Majorana neutrinos to SM particles. 



2.3 Thermal leptogenesis in type-I seesaw models 

In thermal scenario we assume that the out-of-equilibrium decay of the heavy right handed 
Majorana neutrinos to SM leptons and Higgses produce a net lepton asymmetry |HJ. We 
further assume a normal mass hierarchy, Mi -C M 2 < M 3 among the right handed 
Majorana neutrinos Ni,N 2 and N3. In this scenario, while the heavier neutrinos N 2 and 
7V3 decay, the lightest right handed neutrino N% is still in thermal equilibrium. Any L- 
asymmetry thus produced by the decay of N 2 and N3 is erased by the lepton number 
violating interactions mediated by N\. Therefore, it is reasonable to assume that the 
final L-asymmetry is produced by the out-of-equilibrium decays of Ni only. A partial 
L-asymmetry is then transformed to 5-asymmetry via the high temperature sphalerons 
which are in equilibrium at a temperature 10 12 GeV to 10 2 GeV. Below the electroweak 
phase transition (T ew ~ lOOGeV), the sphaleron transitions fall quickly out of thermal 
equilibrium. Therefore, the B-asymmetry produced until at the scale of lOOGeV is the 
final .B-asymmetry of the universe that is observed today. 

2.3.1 Upper bound on CP- asymmetry 

Beyond SM all the processes mediated by right handed neutrino, Nr naturally violate 
lepton number. Out of which the dominant channel is the decays of Nr to SM lepton (£) 
and Higgs (0) through the Yukawa coupling. The required Lagrangian for the coupling is 
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given by 

C Y = I'ijdoSn, + h.c. , (2.6) 

where hij is the Yukawa coupling matrix, and i, j = 1, 2, 3 for three flavors. 

We shall work in a basis in which the Majorana mass matrix of right handed neutrinos 
is diagonal. In this basis the right handed Majorana neutrino is given by Nj = iVjy ±iV^-, 
which satisfies iVJ = ±iVj. The type-I seesaw mechanism then gives the corresponding 
light neutrino mass eigenstates z/ 2 , with masses mi, m 2 , m 3 , respectively; these are 
mixtures of flavor eigenstates v e , v T . 

The decays of the heavy right handed Majorana neutrino can create a non-zero L- 
asymmetry only if their decay violates CP. The CP- asymmetry parameter in the decay 
of Nj is defined as 

3 Y{Nj -> £0) + T(iV, -> £ C C ) ' 1 ' 

In the hierarchical scenario it is reasonable to assume that the decays of N\ is the leading 
process that produces the final L-asymmetry. In this scenario, equation (|2.7jl gives rise to 
the CP- asymmetry, obtained by superposing the tree diagram with the one loop radiative 
correction jH] and self energy correction diagrams to be 




(a) (b) (c) 



Figure 2.1: (a)tree level (b) one loop correction (c) self-energy correction diagrams for 
the decay process of N\. 



1 / m 2 \ 

where g(x) = —3/2y/x in the limit x = (Mj/Mf) ^> 1 (which indicates a hierarchy among 
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right handed neutrinos). In this approximation equation (|2.8|) translates to 



167r[/it/i] n 

3Mi 
167r[/it/ l ] u 



2=2,3 



M 



(2.9) 



Using ()2.3|) and ()2.5|) in equation ()2.9|) we get 

3Mi 



^m^ ,m(h}ULDmUlh " )u 



3Mi 



167ri; 2 [/it/ l ] 11 
We now define EH |2E] 



miJm [(C4/i)n +^milm {U j L h) a 



R = vD; r y 2 UlhD-J /2 , 



(2.10) 



(2.ir 



such that i?i? T = 1. This implies that Im(RR T )u = 0. Then it is straightforward to 
show that 

— Im(Ulh) 2 u = -Im^—ImiUlh) 2 ,. 



Using (|2.12|) in equation (j2.1(Jj) we get 

3Mi x - (mf - m\) 



16TTV*[tfh] U 



rrii 



Im(h 



(2.12) 



(2.13) 



where h = U\h. 



The recent atmospheric neutrino data ^Hl indicates «-> z/ r oscillation with nearly 
maximal mixing(# atm ~ 45°) and a mass-squared-difference 



Ami = K - m l\~ 2 - 6 x W 3 eV 2 . 



(2.14) 

We assume a normal hierarchy among the light neutrino mass eigen states. In this scenario 
the atmospheric neutrino mass is 



V 7 Am 2 atm = m 3 ^ 0.05eV. 



(2.15) 



Using (|2.15|) in equation (J2.13|) we get the upper bound on the CP-asymmetry paramet 
to be 

3Mi 



er 



El < 



16ttv 2 



VAm 2 atm . 



(2.16) 
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2.3.2 Analytical estimation of L-asymmetry and lower bound 
on the mass of lightest right handed neutrino 

In this section, we analytically estimated the L-asymmetry from the out-of-equilibrium 
decays of Aq. We assume that above the mass scale of Aq, all the lepton violating processes 
mediated by Aq are fast enough to bring it in thermal equilibrium. Hence there is no net 
L-asymmetry. Below the mass scale of Aq all these processes fall out of equilibrium. Any 
L-asymmetry thus produced by the decays of Aq is not erased. 
In a comoving volume, the net L-asymmetry is defined by 

Y L = ^ei6, (2.17) 
s 

where s = (2n 2 / '45) g*T 3 is the entropy density at any epoch of temperature T, e\ is the 
CP-violating parameter and 5 is the wash out factor due to the L-violating processes 
mediated by Aq at a temperature T < Mi. A part of the produced L-asymmetry is then 
transferred to L?-asymmetry via the sphaleron transitions and is given by ^1] 

Y B = ^Y L , with P= **F + 4N* 
B p-1 ' y 22A F + 13A/ V ; 

where Np is the number of generation for right handed neutrinos and N<j, is number of 
doublet Higgs. For three generation of right handed neutrinos, p = 28/79. Using (j2.17J) 
and (|2.16|) in equation (J2.18)) we get a bound on the net L?-asymmetry to be 

Y B < ^^i^VA^. (2.19) 
p — 1 s lbirv z 

Recent observations from Wilkinson Microwave Anisotropy Probe {WMAP) show 
that the matter-antimatter asymmetry in the present Universe measured in terms of ^ a 
is 



(—) = (6.118:5) x io- 10 , (2.20) 

V n i J o 

where the subscript presents the matter-antimatter asymmetry that is observed today. 

To compare with the observed value, we now recast the bound on L?-asymmetry (J2.19|) 
in terms of (ng/n 7 ). This is given to be 

(n B /n 7 ) = 7Y B < ^^^V^^*- (2-21) 
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Comparing (J2.2U|) and (J2.21|) we get the constraint on the mass scale of iV\ to be 

M, > 0.84 x l^GeV ( A ( ( ° 05 ^ ) , (2.22) 

1,6.4 x 10-" J \^S) \mGeVJ ' 

where we have used v = 17AGeV, the scale of electroweak phase transition. Other phys- 
ical quantities, the atmospheric neutrino mass and the observed baryon asymmetry, are 
normalized with respect to their observed values. 



2.3.3 Numerical estimation of L-asymmetry 

The analytical estimation in section 2.3.2 shows that in a thermal scenario to create a 
successful L-asymmetry from the decay of Ni it is required that Mi > O(10 9 )GeV . We 
now check the compatibility of this bound on Mi numerically by solving the Boltzmann 
equations El ESI- 

For demonstration purpose, we consider a model [Hlj based on the gauge group 
SU(2) x U(1)y x U(1)y', where Y' is a linear combination of Y and B — L. Since 
B — L is a gauge symmetry, no primordial L asymmetry exists. A net L-asymmetry is 
created dynamically after the B — L gauge symmetry breaking phase transition. A part of 
this asymmetry is then transformed to L?-asymmetry via the non-perturbative sphaleron 
processes. 

In a diagonal basis the right handed Majorana neutrinos acquire masses Mj = fiTjB-L, 
Vb-l being the scale of B — L symmetry breaking phase transition and / being the 
Majorana Yukawa coupling matrix. Above the mass scale of N\ all the interactions 
mediated by N\ are fast enough to keep it in thermal equilibrium. This implies that 
there is no net L-asymmetry. However, as the temperature of thermal plasma falls and 
becomes comparable with the mass scale of Ni, an L-asymmetry is created through the 
CP- violating decays of Ni. However, a part of the created asymmetry is erased by 
the inverse decay and L- violating scatterings mediated by Ni. We study the dynamical 
generation of a net L-asymmetry by solving the relevant Boltzmann equations 



dY 



Nl 



dZ 
dY B _ L 



dZ 
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-(D + S) (Y m - Y e N \) (2.23) 
-tiD (Y m - Y e N \) ~ WY B . L} (2.24) 



where Y a = (n a /s) is the density of any species 'a' in a comoving volume and s is the 
entropy density. Here Z = Mi/T is a dimensionless parameter, where T is related to the 
cosmic time t through the time temperature relation 

t = 0.3g;^. (2.25) 

A derivation of these equations is given in appendix A. The terms D, S and W occurring 
above are now discussed. 

1 . D = accounts for the decay and inverse decay of lightest right handed neutrino 
into lepton and Higgs, Aq <-> 0(0) + 1(1), where 

r * = tL^ m W) (2 - 26) 



r ID = ^T D . (2.27) 
In the equation ()2.26|) the parameter rh\ is defined as 



(m r D m D ) 



*> = (2.28) 

called the effective mass of the light neutrino jTTj. K\ and K 2 are modified Bessel 
functions whose ratio in equation ()2.26|) gives the time dilation factor. At a tem- 
perature above the mass scale of N\ one can check that ~ Tjd- Below its mass 
scale the inverse decays are blocked. So the density of Ni changes significantly due 
to the decays of Aq. 

2. S = accounts for the lepton number violating AL = 1 scatterings. The possible 
reactions Aq+Z(Z) <-> i(t)+q(q) via the exchange of in the s-channel and Ni+t(i) <-> 
/(/) + g(g) through the exchange of in the t-channel, are shown in figure 2.2. The 
total rate of AL = 1 violating scatterings is given as 

T s = 21$ + 4T&, (2.29) 

where 
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(a) 



(b) 



Figure 2.2: AL = 1 lepton number violating scatterings mediated by Standard Model 
Higgs through s or t-channel. 
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Mfrhi 



dx\fxK\{Z *J~x) 



x — 1 + y' 

y' 
l " 



X 



X 



(2.30) 
(2.31) 



128ir 3 v 4 K 2 (Z) j, 

Here we have used y' = (m^/Mf) and the dimensionless quantity x = s'/Mf, with 
s' being the square of center of mass energy. Note that in the above scattering rates 
we have neglected the corrections due to second and third generation right handed 
heavy Majorana neutrinos. 

3. W — constitute the wash out processes which compete with the decay term that 
actually produce the asymmetry. Here Yw receives the contribution from the inverse 
decay (T ID ), AL = 1 scatterings (r}J, rJJ) and AL = 2 scatterings (1^, T l Nt ). The 
AL = 2 scattering processes, £<fi <-> £(f) via the exchange of Ni and ££ <-> 00 mediated 
by Ni in the t-channel, are shown in the figure 2.3. Combining all these processes 
we get the total scattering rate for the wash out processes to be 



1 



rv - t^id + 2 (r l N1 + 2T l Nht ) + 2r^ + -^r^. 



(2.32) 



Nl 



19 



(a) 



(b) 



Figure 2.3: AL = 2, lepton number violating scatterings mediated by lightest right handed 
neutrino. 
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where 



1 



x 



, , with ci = — — 

(a;-l) 2 + ci' M x 



(2.33) 



(2.34) 



(2.35) 



The quantities are thermally averaged reaction rates per particle X. They are related 
to the reaction densities ji as 



n{z) 



n 



cq 



(2.36) 



The reaction densities are obtained from the reduced cross-sections <fi(s'/Mf) as follows- 

M 4 1 



li{Z) 



64vr 4 Z 



dx a i {x)\fxKi{Zy/x) ) 



(2.37) 



<(ml+ml)/M* 

where m a and m& are the masses of the two particles in the initial state. 

Equations (|2.23|) and ()2.24|) have been solved numerically. We assume that far above 
its mass scale the species N\ is in thermal equilibrium. So the initial abundance of Ni 
is determined by its equilibrium distribution. Further at equilibrium, decays or lepton 
number violating scatterings will not produce any asymmetry. Therefore, we assume the 
following initial conditions 



yin _ Y e 1 arid V in — fl 



(2.38) 
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Figure 2.4: Dynamical evolution of neutrino number density in a comoving volume and 
the produced B — L asymmetry at Mi = 10 w GeV and t\ = — 5 x 1CT 8 . 



In fig. 2.4 (a) we have used rh\ = 10~ 6 eV. Far above the mass scale of N% (i.e. Z <§C 1) 
the B — L asymmetry is zero. As the temperature falls (i.e. Z increases) the asymmetry 
builds up and finally it reaches a constant value Yb-l — 2.6 x 10 -10 when the wash out 
effects become negligible. On the other hand, in fig. 2.4 (b) we have used fh\ = 10~ 4 eV 
and we get a smaller asymmetry 1.8 x 1CT 10 . This is because of larger effective neutrino 
mass. Note that the wash out effects not only depend on M% but also rh\. Therefore, for 
a larger rh\ the wash out effects is more and thus we get effectively a smaller asymmetry. 

We now introduce a new parameter m* to be called equilibrium neutrino mass. It is 
defined by jlS] 

m* = (2.39) 
ti 

where 

f d = T7T~2^ M l ( 2 - 4 °) 
is obtained from equation ()2.26|) in the limit Z — > oo and H is the Hubble expansion 
parameter. At any epoch of temperature T, the Hubble expansion parameter is defined 
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as 

H = 1-67^ 1/2 ^. (2.41) 
Using (|2.41j) in equation ()2.39j) we get the equilibrium mass at T ~ Mi to be 

G 1/2 

m„ = 16tt^ 1/2 n -2x 10- 3 eV, (2.42) 
v 1Gp 

where Gn and Gf are Newton and Fermi constants respectively and therefore m* may 
also be called cosmological neutrino mass |43j . 

We now define a dimensionless parameter K — which determines whether the 
species Aq is in thermal equilibrium. For K > 1, the inverse decay processes are fast 
enough to ensure the species Aq to be in equilibrium, irrespective of its initial abundance, 
in the epoch Z — > (T — > oo). In this regime, any pre-existing asymmetry gets erased 
by the rapid inverse decay processes. Therefore, it is called strong wash out regime [2H] • 
In this case the final lepton asymmetry doesn't depend on the initial conditions 1)2.38)1 . 
On the other hand, if K < 1, then the inverse decay processes are suppressed and the 
abundance of Aq is not brought to equilibrium even for Z —>■ 0. In this case, any pre- 
existing asymmetry produced at the B-L symmetry breaking scale continue to be as it is 
until it gets some comparable contribution from the decays of Aq. So this regime is called 
weak wash out regime. In this regime the final lepton asymmetry strongly depends on the 
initial conditions. 

In the following we study the solution of Boltzmann equations (J2.23|) and (|2.24j) by 
taking the zero initial abundance of Aq for two different values of fh\ in the case K < 1. 
For the numerical solution we assume that 

y*J = and Yp_ L = (2.43) 

In the figure 2.5(a) we have used rh\ = 10~ 6 eV, which is three orders of magnitude 
less than the equilibrium mass, m* ~ O(10~ 3 )eV of the light neutrino. In this scenario, 
the lightest right handed neutrino Aq decays being out of equilibrium through out the 
evolution. Hence the erasure any pre-existing L-asymmetry is prevented. This remark 
is relevant to our study in chapter 4. On the other hand, in figure 2.5(b) we have used 
a larger fhi = 10~ 4 eV, which is one order of magnitude less than m*. Therefore, the 
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Figure 2.5: Dynamical evolution of neutrino number density for K < 1 and = along 
with the produced B — L asymmetry at M\ = 10 10 GeV and ei = — 5 x 10 -5 . 



neutrino abundance reaches the equilibrium value at an earlier time than the previous 
case. 

Similar calculations are done for the case K > 1 |2Hj- It is shown that the species 
Ni is brought to equilibrium quickly even if we start with the zero abundance of N%, and 
hence erasing any pre-existing asymmetry, in the epoch Z — > 0. 
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Chapter 3 

Thermal leptogenesis in type-II 
seesaw models and bounds on 
neutrino masses 

3.1 Introduction 

In the type-I seesaw models the upper left 3x3 block of the 6x6 neutrino mass matrix 
is zero; see e.g., section 2.2. This is because of the absence of ul — vl interaction in 
the SM Lagrangian as it violates the lepton number by two units. In contrast to it, in 
type-II seesaw models the presence of an additional scalar triplet A L allows us to add a 
ul — ul interaction to the SM Lagrangian by compensating the two units of B — L charge 
appearing in the ul — ul interaction term. At a low scale the Al acquires a VEV, thus 
providing an additional mass Ml = /(Al), f being the Majorana Yukawa coupling, to 
the light neutrino mass eigenstate through the diagonalization of the 6x6 neutrino mass 
matrix. As a result the light neutrino mass matrix takes the form m v = M L — m 2 D /M R . 
The two terms are called type-II and type-I respectively. The class of models in which 
both type-I and type-II terms occurring in m v are called type-II seesaw models. 
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3.2 Type-II seesaw mechanism and neutrino masses 

In the minimal scenario, to achieve the light neutrino masses via the type-II seesaw 
mechanism, a scalar triplet and a right handed Majorana neutrino per family are 
added to the SM. Thus the Lagrangian of this model reads 

£ = C-SM + C-new (3.1) 

where Csm is the SM Lagrangian and C new is the additional Lagrangian that contains 
the new interaction involving the right handed neutrinos and the triplet A^. The relevant 
terms of the Lagrangian are given to be 

-Cnew ^ -MRiN^CNRi + MlTrAlAL+hiji^NRj 

+ fijfCi^Aij - /j(f) T iT 2 A L (j) + H.C. (3.2) 

where i T = [y Li -, &ia) and <p T = (0°, <fi~) and 

T, 5+ s++ 



(3.3) 

In the presence of these interactions, the neutral component of the triplet acquires a VEV, 



, 0\ 

(Al> =| , (3.4) 

v L 0) 

at a scale much below the electroweak symmetry breaking phase transition. Due to this 
VEV there are now in general two sources of light neutrino masses 

m v = fv L -v 2 hM R 1 h T 

= 171^ + ml (3.5) 

Note that m j that contribute to the neutrino mass matrix in the present case was absent 
in type-I models. 

We can diagonalize the light neutrino mass matrix m u , through lepton flavor mixing 
matrix Ul [HI]- This gives us three light Majorana neutrinos of masses 

Ulm u U L = dia(m 1 ,m 2 ,m 3 ) = D m , (3.6) 

where the masses mi,m 2 and m 3 can be chosen to be real and positive. 
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3.3 Thermal leptogenesis in type-II seesaw models 



In the type-II seesaw models the following decay modes: 



A 




f + f 



and 



N 




1 + <J> 



violate lepton number by two units and hence produce the lepton asymmetry. In the 
above equations I and <fi are SM lepton and Higgs. 

In what follows we assume a normal mass hierarchy in the heavy Majorana neutrino 
sector. We also assume that the quartic self coupling of the SM Higgs, which is expected 
to be of order unity, is much larger than the Majorana Yukawa coupling of lightest right 
handed heavy Majorana neutrino N\. In this case while the heavier right handed neutri- 
nos, N 2 and iV 3 and the triplet A^ decay, the lightest of heavy Majorana neutrinos is still 
in thermal equilibrium. Any asymmetry thus produced by the decay of N 2 , iV 3 and Al 
will be erased by the lepton number violating interactions mediated by Ni. Therefore, it 
is reasonable to assume that the final lepton asymmetry is given only by the CP- violating 
decays of Ni to the SM fields £ and 4>. 

3.3.1 Upper bound on CP-asymmetry 

In comparison to the CP-asymmetry (|2.8|) in type-I models there is an additional contri- 
bution EE] in type-II seesaw models due to the one loop radiative correction through 
the virtual triplet A^ in the decays of lightest right handed Majorana neutrino. We as- 
sume that the masses of Al, N 2 and N 3 are much heavier than the the mass scale of Ni. 
In this scenario the total CP-asymmetry is given by 



where the contribution to e[ comes from the interference of tree level, self-energy correction 
and the one loop radiative correction diagrams involving the heavier Majorana neutrinos 



ei = e[ + e 



ji 

1 j 



(3.7) 
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N 2 and N 3 . This contribution is the same as in type-I models [2ZH2B] an d is given by 
On the other hand the contribution to e{ 7 in equation ()3.7j) comes from the interference 



4» 

y\ — » 




Figure 3.1: One loop radiative correction through the virtual triplet A^. 



of tree level diagram and the one loop radiative correction diagram involving the triplet 
Ai as shown in fig. 3.1. It is given by |31| 167] 

Substituting ()3.8|) and (j3.9|) in equation (|3.7|) and using (|3.5|) . we get the total CP- 
asymmetry 

3Mi Im(tfm v h*) u in . 



Using ()3.6|) in equation ()3.10|) we get 



3Ml Im(tfUlD m W L h*)u 

Cl 167TU 2 (/lt/l) u 

3M X Ej^MgjMl! rqm 
w EJ(^)a| 2 • 1 • j 

With an assumption of normal mass hierarchy for the light Majorana neutrinos the upper 
bound on CP-asymmetry (|3.11|) can be given by 

3Mi 

ex < jTr^rns. (3.12) 

L07TV 

Note that the above upper bound ()3.12j) for ei as a function of Mi and rris was first 
obtained for the case of type-I seesaw models [23 ■ However, the same relation holds in 
the case of type-II seesaw models also [3T] independent of the relative magnitudes of ml 
and mj/. 
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3.3.2 Estimation of L-asymmetry and lower bound on the mass 
of lightest right handed neutrino 

Assuming Mi <C M 2 , M 3 , Ma, the final lepton asymmetry is given by the out of equilib- 
rium decays of the lightest right handed Majorana neutrino Aq. A part of this asymmetry 
is then transformed to i?-asymmetry by the thermally equilibrated sphaleron processes 
which violate B + L quantum number of SM fermions. In a comoving volume a net 
5-asymmetry can be written as 

Tin 

Y B = — = 0.5561^5, (3.13) 
s 

where the factor 0.55 in front |Hj is the fraction of L-asymmetry that is converted to 
5-asymmetry. Here Y^ 1 is the density of Aq in a comoving volume which is given by 
Yzvj — n Ni_/s, s being the entropy density of the Universe at any epoch of temperature T 
and S is the wash out factor. 

We now recast equation (j3.13|) in terms of a measurable quantity (ne/w-y) which is 
given by 

(j^J = 7(Y B ) = 3.85(61*^). (3.14) 
Substituting equation (|3.12j) in (|3.14j) we get a bound on the baryon asymmetry to be 

-) < 3-85 (^) m 3 Y Nl 5. (3.15) 

Using ()2.15j) in equation (j3.15j) and comparing with the observed baryon asymmetry ()2.20|) 
we get a bound on the mass of Aq to be 

Ml > 0.84 x W °GeV ) . (3.16) 

\Y Nl SJ \ m 3 J 

This bound was obtained in type-I seesaw models. However, in this case we revisit the 
same bound on the mass of Aq irrespective of any assumption regarding the magnitude 
of type-I and type-II terms in the neutrino mass matrix f)3.5j) . 



28 



3.4 Spontaneous CP- violation and leptogenesis in Left- 
Right symmetric models 

In section 3.2 we demonstrated the type-II seesaw mechanism in a minimal scenario by 
adding a right handed Majorana neutrino per generation and a heavy triplet Al to the 
SM. However, the light neutrino masses via type-II seesaw mechanism can be obtained 
naturally in Left-Right or SO (10) models. In the following we consider the low energy 
left-right symmetric model in which we assume the case of spontaneous CP-violation 
(SCPV). In this scenario we derive an upper bound on the CP-asymmetry. Moreover, 
we discuss the bounds on neutrino masses from the leptogenesis constraint. 

3.4.1 Left-Right symmetric model and SCPV 

In the low energy left-right symmetric model the right handed charged lepton of each 
family, which was a singlet under the SM gauge group SU(2) L <g> U(l)y, gets a new 
partner v R . These two form a doublet under the SU(2) R of the left-right symmetric gauge 
group SU(2) L <g> SU(2) R <g> U(1) B -l- Similarly, in the quark sector, the right handed up 
and down quarks of each family, which were singlets under SM gauge group, combine to 
form a doublet under SU(2) R . 

The Higgs sector of the model is dictated by two triplets A L and A R and a bidoublet 
$, which contains two copies of SM Higgs. Under SU(2) L ® SU(2) R ®U(1) B -l the field 
content and the quantum numbers of the Higgs fields are given as 



To achieve the correct phenomenology, the various Higgs multiplets in the model 




(3.19) 



(3.17) 



(3.18) 
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should have the following VEVs, 




(3.20) 



(3.21) 



(3.22) 



Q = Tl + Tl + -{B-L). 



(3.23) 



In the above vl, vr, k\ and &2 are real parameters and the electroweak symmetry breaking 
scale v = 174 GeV is given by v 2 = k\ + k\. Further we require that vl <C v <C Vr. The 
requirement of the spontaneous breakdown of parity gives rise to 



where 7 is parameter which is a function of the quartic couplings in the Higgs potential. 

The minimisation of the most general Higgs potential involving Al,Ar and $ was 
studied in refs. [HI]. The relations between the various couplings, for which the above set 
of VEVs are generated, were derived. In this scenario, the gauge symmetry SU(2)l x 
SU{2)r x U(1)b-l is broken to U(l) em in a single step. Thus the CP-violating phases 
come into existence at the same scale where the left-right symmetry is broken. Since 
v Vr, the SM symmetry is present as an approximate symmetry at the scale where 
symmetry breaking occurs. 

The fermions get their masses via Yukawa couplings. The Lagrangian for one genera- 
tion of quarks and leptons is 




(3.24) 



yuk 



h q q L <&q R + g q qL$qn + h t £ L ^l R + gi£ L $l R 




(3.25) 
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where q and t are quark and lepton doublets, $ = r 2 $*r 2 and C is the Dirac charge 
conjugation matrix. Further the Majorana Yukawa coupling / is the same for both left 
and right handed neutrinos to maintain the discrete L «-> R symmetry. 

Substituting the complex VEVs flSHD , (EQ) and (j3~22|) in fl£Z5D we obtain fermion 
mass terms to be 

£ ma ss = (h q k x e* a + g q k 2e lP )u L u R + (h q k 2 e tf3 + g q k x e %ol )d L d R 
+ (h l k 1 e ia + gik 2 e il3 )u L u R + (hihe^ + gik ie ia )e L e R 
+f(v T L Cv L e te ^ L + ulCv R e i0 *v R ) + H.C. (3.26) 

Generalizing the above equation ()3.2fij) for three generation of matter fields we get the up 
and down quark mass matrices to be 

(M u ) i3 = (h q hhe ta + (g q )ijk 2 e ip and (M d ) t] = (hJufoe* + (g^jhe™. (3.27) 

We assume j^HEHl k\jk 2 ~ m t /mb. In the seesaw mechanism, the Dirac mass matrix of 
the neutrinos is assumed to be similar to the mass matrix of the charged leptons. For 
k 2 -C and further assuming hi ~ gi in (|3.2fjjl . the Dirac mass matrix of the neutrinos 
to a good approximation becomes hikie lOL . Thus neglecting k 2 terms, the masses of three 
generations of neutrinos are given by 

Cv-mass = HAe^ChihvR, + hiy^vlfv^ + v R e ie R v T R C v R] ) + H.C. (3.28) 

The Majorana mass matrix for the right handed neutrinos can be diagonalized by making 
the following orthogonal transformation on v R 

N R = T R v R . (3.29) 

In this basis, we have 

T R fO R = f dm , (3.30) 

h = hO R . (3.31) 

In the transformed basis we get the mass matrix for the neutrinos 

fv L e if>L k x e ia h 
he ia h T f dia v R e i6 « 
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Diagonalizing the neutrino mass matrix into 3x3 blocks we get the light neutrino mass 
matrix to be 

m v = fv L e»> - ^{hf- dl lh T y^-^ (3.33) 



Notice that the Lagrangian ()3.25|) is invariant under the following unitary transforma- 
tions of the fermion and Higgs fields, 

if) L — > U L ip L and i) R — > U R ^ R , (3.34) 
$ — ► U L $U R and $ — ► U L $U R (3.35) 
A L -^U L A L U ] L and A R — ► UrArU r , (3.36) 

where ipL,R is a doublet of quark or lepton fields. The invariance under Ui is the result 
of the remnant global U(l) symmetry which remains after the breaking of the gauge 
symmetry SU (2)l and similarly for Ur. The matrices Ui and U R can be parametrized as 

( e iiL o \ fe^ R , 

U L =\ and [/^ . (3.37) 

I o e ~^ L I V e~^ R 

By redefining the phases of the fermion fields we can rotate away two of the phase degrees 

of freedom from the scalar sector of the theory. Thus only two of the four phases of Higgs 

VEVs have phenomenological consequences. Under these unitary transformations, the 

VEVs (BID , (I3~2T| and lCP2|) become 

(A R )= I ° °], (3.38) 

( kie i( - a+ ^ L ~^ R) 

A; 9 e i(/3 - 7i+7 « ) 



and 



(3.39) 



\ 

(AO = | • (3.40) 

Uzel (fe-2 7i ) o y 

We choose 7^ = 0r/2 so that the masses of the right handed neutrinos are real. The light 
neutrino mass matrix ()3.33|) then becomes 

„2 



V R 

^ + mi 



mi 1 + mi (3.42) 
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Conventionally, in equation (|3.41|) . ^ L was chosen to be —a-\-9n/2 [HUES]- This makes 
ml real leaving the imaginary part purely in m}J . We call this type- II phase choice. The 
light neutrino mass matrix, with this phase choice, is 

m v = fv L e^ - f(hf^ a h T ), (3.43) 

where 9' L = {9 L — 9 R + 2a). On the other hand, by choosing jl = 9l/2 in equation ()3.4H) 
m^ can be made real, with the phase occurring purely in m} v . We call this type-I phase 
choice. Consequently the light neutrino mass matrix (|3.41j) becomes 

m v = fv L - —e i0 R(hf^lh T ) (3.44) 

where 9' R = (9l — 9r + 2a). The CP- violating parameter e\ which gives rise to the lepton 
asymmetry is independent of the phase choice. However, the theoretical upper bound on 
ei is not a physical parameter of the theory and can depend on the choice of phases as we 
see in the next section. In numerical calculations, we take into account the consistency 
of the bounds coming from the different phase choices. 

Using ()3.6|) we can diagonalize the light neutrino mass matrix m v . This gives us three 
eigenvalues, mi, 7712 and 777,3 which are chosen to be real. 

3.4.2 Upper bound on CP- asymmetry in Left-Right symmetric 
models with SCPV 

Following the same convention in section 3.3.1 we can write the total CP-asymmetry in 
Left-Right symmetric model as 

3Mx Eij Im [hUml + ml^hj!] 

61 = Tr — 2 7WT\ ' 3 - 45 

lbnv z [h 1 n)n 

From equation (|3.45|) . we see that the physical observable ei is not affected by the choice 
of phases. In the following, we use bound on ei from the observed baryon asymmetry to 
obtain bounds on right-handed neutrino masses for the two different phase choices. 

A. The type-II choice of phases 

In this choice of phases the type-I mass term is real. The only source of CP-violation in 
the light neutrino mass matrix m u lies in the type-II mass term. Thus in this case e{ = 
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because of both h and m T u are real. The total CP-asymmetry in this choice of phases is 
therefore given by 



3M lVL (h T fh) 



11 



lQirv 2 (h T h) u 
Using ()3.3()j) and (j3.31j) in equation ()3.46j) we get 



Jm(e^). (3.46) 



f _ZM 1 v L Y. l h{Q T R h)l 
IGttv 2 22il°Rh)a 

where fa = (Mj/Vr). Up to a first order approximation it is reasonable to assume that 
^\ /j pa 1. In this approximation the maximum value of the CP-asymmetry ()3.47|) is 
given by EH 123 ED 



£l,max = TZZ^t ' (3.48) 



3M lVL 

167TV 

Thus, for type-II phase choice, a bound on e\ leads to a bound on M\. 
B. The type-I choice of phases 

In the type-I choice of phases the type-II mass term is real. Hence the CP-violation comes 
through the type-I mass term only. The total CP-asymmetry in this case is therefore given 
by 



ei = A 



3M l k 2 (h T hf-lh T h)n T _^_ ld 



16ttv 2 vr (h T h)n 
Let us consider the type-I term of the light neutrino mass matrix 



Im(e- iy K). (3.49) 



ml = m v — ml 

2 



~^hf^h T , (3.50) 

U R 

where v' R = vrc^'r. We can find a diagonalizing matrix U = OU p h aS e for ml such that 

U T mlU = —D mi = —dia(m,i r ,mi 2 ,mi^) (3.51) 

where mj 1 ,mj 2 and mi 3 are made real by choosing U p hase = e 1 ®'^ 2 . Therefore, from 
equation (|3.51|) we have 

D mi = k lo T (hf^ a h T ) O. (3.52) 
Vr 
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Using (|3.52|) in equation ()3.49|) the CP-asymmetry e\ can be rewritten as 

3Mi 



1671V 2 E t [(h T 0) ll (0 T h) a ] 
3Mx £^(^)a /m(e -< 



Im(e 



-ie' 



(3.53) 



In the above equation 1)3.53)1 the maximum value of CP- asymmetry is thus given by |2"?|l2"£] 



e i,' 



3Mi 



E 



(3.54) 



In the equation ()3.54j) mjs are the eigenvalues of the matrix ml and are not the physical 
light neutrino masses. It is desirable to express the e^ max in terms of physical parameters. 
In order to calculate the m/S we will assume a hierarchical texture of Majorana coupling 



Mi 

Jdia 

Vr 



\ 



(3.55) 



1 

OCA 

\0 a B J 

where 1 < = (M 2 /Mi) <C as = (M3/M1). We identify the neutrino Dirac Yukawa 
coupling /i with that of charged leptons [TH] . We assume h to be of Fritzsch type [70] 



(m T /v) 
1.054618 



A) a 0^ 



a b 
6c 



(3.56) 



We make this assumption because Fritzsch mass matrices are well motivated phenomeno- 
logically. By choosing the values of a, b and c suitably one can get the hierarchy for 
charged leptons and quarks. In particular [7UJ 



a = 0.004, b = 0.24 and c = 1 



(3.57) 



can give the mass hierarchy of charged leptons. For this set of values the mass matrix 
h is normalized with respect to the r-lepton mass. The set of values of a, b and c are 
roughly in geometric progression. They can be expressed in terms of the electro-weak 
gauge coupling a w = = si " 2e . In particular a = 2.9a 2 , b = 6.5a,,, and c = 1. Here 
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onwards we will use these set of values for the parameters of h. Using equation (J3.55)) and 
(I3.56J) in equation ()3.52j) . we now get 

„,2 



D 



ml 



Vr 



dia 



(m T /vy 



Mi (1.054618) 



/r 









.4 








B 



(3.58) 



where the eigenvalues A and B are functions of a a and o>b and their sum is given by 

A + B= 1 -\a 2 + —(a 2 + b 2 ) + — 
2 [ a a O-b 

Using equation ()3.58|) we can write the maximum value of CP-asymmetry ()3.54|) 

3Mi 



(3.59) 



16ttv 2 
3 (m T /v) 2 



m I2 + m I3 ) 

(A + B) 



(3.60) 



16tt (1.054618) 2 

Thus we see that, in type-I choice of phases, the leptogenesis parameter ei constrains the 
hierarchy parameters a a and as- In the following two sections, we will obtain numerical 
bounds on a a and «b ina manner consistent with the bound M\ coming from the type-II 
phase choice. 



3.4.3 Estimation of L-asymmetry and bound on neutrino masses 

A net B — L asymmetry is generated when left-right symmetry breaks. A partial B — L 
asymmetry is then gets converted to 5-asymmetry by the high temperature sphaleron 
transitions. However these sphaleron fields conserve B — L. Therefore, the produced 
B — L asymmetry will not be washed out, rather they will keep on changing it to B- 
asymmetry. Thus in a comoving volume a net 5-asymmetry is given by 

Y B = ^ = ^ £l Y m S, (3.61) 

where the factor || in front [H] is the fraction of B — L asymmetry that gets converted 
to fi-asymmetry. Here e\ is given by equation (|3.60p . The other symbols involved in 
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equation (|3.61|) carry the usual meaning; see, e.g. section 3.3.2. However, the observed 
baryon asymmetry of the Universe is measured in terms of ris/n^. Therefore, we rewrite 
equation (j3.61|) as 

^) = 7{Y B ) = 2.4B(e 1 Y N S). (3.62) 

Substituting the type-II phase choice relation ()3.48|) in ()3.62|) and comparing with the 
observed value ()2.20|) of the baryon asymmetry we get the bound on the mass of lightest 
right handed neutrino to be 

M x > 1.25 x 10 8 GeV ( ( ^^-) . (3.63) 

\Y Nl 5J \ v L J 

On the other hand, substitution of €i tTnax from the type-I phase choice (J3.6U)) in equation 
()3.62|) and then comparison with the observed value ()2.20|) gives the constraint 

A + B > 3.46 x UTWrf (j^) (^J , ,3,4, 

where the physical quantities are normalized with respect to their observed values. The 
above equation, for the values of a, b and c from (|3.57|) . gives only one constraint on 
the two hierarchy parameters a a and a B - We will determine the individual parameters 
a a and a B by demanding that their values should reproduce the low energy neutrino 
parameters correctly, while satisfying the inequalities Mi > O(10 8 ) GeV and «b > >> 
1. Individual bounds on a a and a B can also be obtained if we assume that the a a term 
and the o>b term in the sum A + B from equation (|3.59j) are roughly equal. We then get 

a A = (M 2 /Mi) < 17 and a B = (M 3 /M 1 ) < 289. (3.65) 

3.4.4 Checking the Consistency of /-matrix eigenvalues 

The solar and atmospheric evidences of neutrino oscillations are nicely accommodated in 
the minimal framework of the three-neutrino mixing, in which the three neutrino flavors 
z/ e , h>n, v T are unitary linear combinations of three neutrino mass eigenstates u±, i>2, v% with 
masses mi, mz, rn 3 respectively. The mixing among these three neutrinos determines the 
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structure of the lepton mixing matrix [61 which can be parameterized as 



/ 



U, 



cic 3 sic 3 s 3 e'° 

' s l c 2 ~ c l s 2 s 3 e ^ c l c 2 — s l s 2 s 3 e ^ S 2 C 3 

id 



\ 



dia(l,e ia ,e^ +s) ), 



(3.66) 



\ S \S2 — C1C2S3 — C1S2 — SiC2S^e C2C3 J 

where Cj and Sj stands for cos #7 and sin9j. The two physical phases a and (3 associated 
with the Majorana character of neutrinos are not relevant for neutrino oscillations [71] 
and will be set to zero here onwards. While the Majorana phases can be investigated 
in neutrinoless double beta decay experiments [12], the CKM-phase 5 G [— 7r, 7r] can be 
investigated in long base line neutrino oscillation experiments. For simplicity we set it 
to zero, since we are interested only in the magnitudes of elements of Ul- The best fit 
values of the neutrino masses and mixings from a global three neutrino flavors oscillation 
analysis are [13] 



9 1 = 9 Q ^3A°, 9 2 = 9, 



at 111 



45°, fl 3 < 10 c 



(3.67) 



and 



A 2 2 2 

Am = m 2 — m 1 



Am 



at in 



Tfln 



m 



7.1 x 1(T 5 eV 2 
2.6 x 10" 3 eV 2 . 



Using equation (J3.33j) we rewrite the /-matrix 



(m u /eV) 



;i.054165) 2 (Mi/GeV) 



a 2 + 



ab 
be 

Q-B 



ab be b 2 1 (?_ 

oca a B oa a B 



(3.68) 



(3.69) 



where the neutrino mass matrix m u is given by equation ()3.6J) . The constrained eigenvalues 
a a and are given by equation (|3.65|) . 

In the following, we choose M\ to be larger than the bound given by type-II phase 
choice ()3.63|) and m\ such that m\ « A so i. For such mi and Mi, we choose suitable a a 
and as that are compatible with the low energy neutrino oscillation data. In particular 



here we choose m 1 = 1.0 x lO~ 3 eV, M 1 = 1.0 x 10 8 GeV, a A = 17, a B = 170 and 9 3 



6°. 
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Then we get 



dia 



2.16 x lQ~ 3 eV 



1 
17.3 
169.7 



(3.70) 



Thus, for the above values of mi and Mi, the assumed hierarchy of right-handed neutrino 
masses is consistent with global low energy neutrino data. Comparing equation (J3.7(Jj) 
with (|3.55|) we get 

M- o 1 fi v 1 r\- 3 ^v 

-. (3.71) 



Mi 



vl 



This implies that Vr = O(10 10 ) GeV for vl = 0.1 eV. These values of Vl and vr are 
compatible with genuine seesaw vlVr = 'jv 2 for a small value of 7 ~ O(10 -4 ) 74J. On 
the other hand, if we choose the parameters ni\ = 1.0 x 10~ 3 eV, Mi = 1.0 x 10 9 GeV, 
a a = 17, as = 65 and 6*3 = 6° we get 

/ 1 n n \ 



dia 



1.6 x 10- 3 eV 



1 
16.76 
64.68 



(3.72) 



Once again we have consistency between the assumed hierarchy of right-handed neutrino 
masses and global low energy neutrino data. Again comparing equation (|3.72j) with (J3.55|) 

we get 

M, 1 P, v 1 n-3^T/ 

-. (3.73) 



Mi 



Thus for vl = 0.1 eV one can get vr = O(10 n GeV). Again these values are compatible 
with seesaw for 7 ~ O(10 -3 ). 

Here we demonstrated the consistency of our choice of the matrix / with neutrino 
data for two different choices of a a and olb- For other choices of these parameters, to be 
consistent with 1 << a. a « o>b, one can choose appropriate values of m\ < 10~ 3 eV 
and Mi > 10 8 GeV in equation ()3.69|) which will reproduce the correct eigenvalues of the 
matrix /. 
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Chapter 4 

Gauged B — L symmetry and upper 
bounds on neutrino masses 

4.1 Introduction 

It has long been recognized that the existence of heavy Majorana neutrinos has important 
consequences for the baryon asymmetry of the Universe |Hj. With the discovery of the 
neutrino masses and mixings, it becomes clear that only B — L can be considered to 
be a conserved global symmetry of the SM and not the individual quantum numbers 
B — L e , B — and B — L T . Combined with the fact that the classical B + L symmetry 
is anomalous O EH IHj it becomes important to analyse the consequences of any B — 
L violating interactions because the two effects combined can result in the unwelcome 
conclusion of the net baryon number of the Universe being zero. 

At present two broad possibilities exist as viable explanations of the observed baryon 
asymmetry of the Universe. One is the baryogenesis through leptogenesis This has 
been analysed extensively in El EB1 and has provided very robust conclusions for 
neutrino physics. Its typical scale of operation has to be high, at least an intermediate 
scale of 10 9 GeV. This has to do with the intrinsic competition needed between the lepton 
number violating processes and the expansion scale of the Universe. While the mechanism 
does not prefer any particular unification scheme, it has the virtue of working practically 
unchanged upon inclusion of supersymmetry [75] . 
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The alternative to this is provided by mechanisms which work at the TeV scale (35] and 
may rely on the new particle content implied in supersymmetric extensions of the SM. It 
is worth investigating other possibilities jlHj, whether or not supersymmetry is essential to 
the mechanism. The starting point is the observation |4*4~] 143] that the heavy neutrinos 
participate in the erasure of any pre-existing asymmetry through scattering as well as 
decay and inverse decay processes. Estimates using general behavior of the thermal rates 
lead to a conclusion that there is an upper bound on the temperature Tb_l at which B — L 
asymmetry could have been created. This bound is T B _ L <, 10 13 GeVx (leV/m^ 2 , where 
m v is the typical light neutrino mass. This bound is too weak to be of accelerator physics 
interest. We extend this analysis by numerical solution of the Boltzmann equations and 
obtain regions of viability in the parameter space spanned by rhi-Mi, where rhi is the 
effective light neutrino mass parameter as defined in eq. ()2.28j) and M% is the mass of 
the lightest of the heavy Majorana neutrinos. We find that our results are in consonance 
with where it was argued that scattering processes provide a weaker constraint than 
the decay processes. If the scatterings become the main source of erasure of the primordial 
asymmetry then the constraint can be interpreted to imply Tb-l < M\. Further, this 
temperature can be as low as TeV range with fh\ within the range expected from neutrino 
observations. This is compatible with see-saw mechanism if the "pivot" mass scale is a 
factor of 10 2 smaller than that of the charged leptons. 

Here we assume that a lepton asymmetry is produced when the B — L gauge symmetry 
is broken without referring to any specific unification scheme. However, in [7.0] it was 
shown that the Left-Right symmetric model [77] presents just such a possibility. In this 
model B — L appears as a gauged symmetry in a natural way. The phase transition 
is rendered first order so long as there is an approximate discrete symmetry L <-> R, 
independent of details of other parameters. Spontaneously generated CP phases then 
allow creation of lepton asymmetry. We check this scenario here against our numerical 
results and in the light of the discussion above. 
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4.2 Erasure constraints: An analytical estimation 



The presence of several heavy Majorana neutrino species (Ni) gives rise to processes 
depleting the existing lepton asymmetry in two ways. They are (i) scattering processes 
(S) among the SM fermions and (ii) Decay (D) and inverse decays (ID) of the heavy 
neutrinos. We assume a normal mass hierarchy among the right handed neutrinos such 
that only the lightest of the right handed neutrinos (JVi) makes a significant contribution 
to the above mentioned processes. At first we use a simpler picture, though the numerical 
results to follow are based on the complete formalism. The dominant contributions to the 
two types of processes are governed by the temperature dependent rates 

1 and r s ~— - — — — 5 T , (4.1) 



16tt(4T 2 + M 1 2 ) 1 /2 ' ~ 6 13tt 3 (9T 2 + M 2 ) ' 

where h is typical Dirac Yukawa coupling of the neutrino. 

Let us first consider the case Mi > Tb-l- For T < Tb-l, the N% states are not 
populated, nor is there sufficient free energy to create them, rendering the D-ID processes 
unimportant. We work in the scenario where the sphalerons are in equilibrium, maintain- 
ing rough equality of B and L numbers. As the B — L continues to be diluted we estimate 
the net baryon asymmetry produced as [TH] 

Tb ~ l T s dT 




ht)' < 4 ' 2 > 

where ts-L is the time of the (B — L)-breaking phase transition, H is the Hubble param- 
eter, and t EW and T EW corresponds to the electroweak scale after which the sphalerons 
are ineffective. Evaluating the integral gives an estimate for the exponent as 

3\/lO h4 M Pl T B -L 
13vr 4 In 10^ M 2 

The same result upto a numerical factor is obtained in [TH], the suppression factor oj^ 
therein. Eq. (J4.3j) can be solved for the Yukawa coupling h which gives the Dirac mass 
term for the neutrino 

h* < 3200 d J - M \ f ) (4.4) 



T B - L M Plj 

where we have taken g* = 110 for definiteness and d B here stands for total depletion 
including from subdominant channels. This can be further transformed into an upper 
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limit on the light neutrino masses using the canonical seesaw relation. The constraint 
(I4.3J1 can then be recast as 

ftl < l^teY*. (4.5) 



VT B - L M Pl \W 

This bound is useful for the case of large suppression. Consider d B = 10. If we seek 
Tb-l ~ ITeV and Mi ~ lOTeV, this bound is saturated for m v ~ 50keV, corresponding 
to h ~ m T /v. This bound is academic in view of the WMAP bound ^2m Ui w 0.69el/ [HB] 
. On the other hand, for the phenomenologically interesting case m v « lCT 2 eV, with 
ft. w 10~ 5 ~ m e /v and with Mi and Tb-l as above, eq. (|4.5|) can be read to imply that in 
fact d B is vanishingly small. This in turn demands, in view of the low scale we are seeking, 
a non-thermal mechanism for producing lepton asymmetry naturally in the range 1CT 10 . 
Such a mechanism is discussed in sec. 14.41 

In the opposite regime M\ < T B _ L , both of the above types of processes could freely 
occur. The condition that complete erasure is prevented requires that the above processes 
are slower than the expansion scale of the Universe for all T > M\. It turns out to be 
sufficient [13] to require V D < H which also ensures that T s < H . This leads to the 
requirement 

£1/2 

mi < to* = 16ngl /2 N ~ 2 x i(T 3 e1/ (4.6) 
y2Gp 



where the parameter to* jUj contains only universal couplings and <?*, and may be called 
the cosmological neutrino mass. 

The constraint of equation ()4.6|) is compatible with models of neutrino mass if we 
identify the neutrino Dirac mass matrix mo as that of charged leptons. For a texture of 
m D to be Fritzsch type (|3.56j) . equation (|2.28j) gives 

toi ~ 2 x 10- 3 eV ( ^-^-) ■ (4.7) 



Mi J 

Thus with this texture of masses, eq. ()4.6|) is satisfied for Mi ^ O(10 7 )GeV . If we seek 
Mi mass within the TeV range, this formula suggests that the texture for the neutrinos 
should have the Dirac mass scale smaller by 1CT 3 relative to the charged leptons. 

The bound ()4.6|) is meant to ensure that depletion effects remain unimportant and is 
rather strong. A more detailed estimate of the permitted values of toi and Mi is obtained 
by solving the relevant Boltzmann equations in a scenario T^~l > Mi and < H. 
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4.3 Numerical solutions of Boltzmann equations 

The relevant Boltzmann equations (HI El EHl EH] for our purpose are 

dY N1 



dZ 
dY B -L 



(D + S) (Y m - Y e N \) (4.8) 
WY B - L ., (4.9) 



dZ 

The terminologies used in the above equations are explicitly elaborated in chapter 1. We 
will summarize later by considering their importance from the view of present context. 
The thermal corrections to the above processes as well as the processes involving the 
gauge bosons may have significance for final L-asymmetry [59j. However their importance 
is under debate [79: • Therefore, we limit ourselves to the same formalism as in El I2H] • 
In equation (|4.8|) D = To/ HZ, where determines the decay rate of N±, S = 
Fs/HZ, where T$ determines the rate of Al = 1 lepton violating scatterings. In equation 
(J4.9|) W = Tyy/HZ, where Y\y incorporates the rate of depletion of the B-L number 
involving the lepton violating processes with Al = 1, Al = 2 as well as inverse decays 
creating Aq. The various T's are related to the scattering densities [D] 7s as given by 
equation (J2.36j) . As the Universe expands these T's compete with the Hubble expansion 
parameter (H). Therefore, for the AL = 1 lepton number violating processes in a comoving 
volume, we have 



Id \ I l<t>,s \ I l^t 
sE{M x ) J ' \sH(Mi) J ' I sfr(Afi) 



oc fcimi. (4-10) 



On the other hand the dependence of the 7's in Al = 2 lepton number violating processes 
on mi and Mi are given by 

7tvi \ / 7f l m,t \ , - 2 



cx k 2 m{Mi. (4.11) 



^sH(Mi) J ' \sH(Mij 
Finally there are also lepton conserving processes where the dependence is given by 



oc hMi 1 . (4.12) 



^sH(Mi) 

In the above equations (|4.10j) . (|4.11|) . (|4.12jh hi, i = 1,2,3 are dimensionful constants 
determined from other parameters. Since the lepton conserving processes are inversely 
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proportional to the mass scale of Aq, they rapidly bring the species Aq into thermal 
equilibrium for Mi < T. (|4.11|) are negligible because of their linear dependence on Mi. 
This is the regime in which we are while solving the Boltzmann equations in the following. 
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Figure 4.1: The evolution of B-L asymmetry for different values of Mi shown against 
Z{= Mi/T) for rhi = 10" 4 eV and r] raw = 2.0 x 1(T 10 

The equations (|4.8|) and (J4.9)) are solved numerically. The initial B — L asymmetry is 
the net raw asymmetry produced during the B-L symmetry breaking phase transition by 
any thermal or non-thermal process. As such we impose the following initial conditions 

Y in — V eq anrl V in — n raw 
1 Nl — 1 Nl dna 1 B-L — 'IB-L 

At temperature T > M 1; wash out effects involving N\ are kept under check due to 
the m\ dependence in (|4,11|) for small values of fh\. As a result a given raw asymmetry 
suffers limited erasure. As the temperature falls below the mass scale of Mi the wash out 
processes become negligible leaving behind a final lepton asymmetry. Fig. 14.11 shows the 
result of solving the Boltzmann equations for different values of Mi. 

If we demand that the initial raw asymmetry is of the order of ax 1CT 10 , with a ~ 0(1), 
then in order to preserve the final asymmetry of the same order as the initial one it is 
necessary that the neutrino mass parameter rhi should be less than 10~ 3 e\A This can be 
seen from fig. 14.21 For rhi — 10~ 3 e\^ we can not find any value of Mi to preserve the 
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Figure 4.2: The allowed values of Mi against the required final asymmetry is shown for 
7] raw = 2.0 x 1(T 10 

final asymmetry, a x 10~ 10 in the allowed region. This is because of the large wash out 
effects as inferred from the equation (|4.11jl . However, for mj = 10 eV we get a lowest 
threshold on the lightest right handed neutrino of the order 10 12 GeV. For any value of 
Mi < 10 12 GeV the final asymmetry lies in the allowed region. This bound increases by 
two order of magnitude for further one order suppression of the neutrino parameter rhi. 
The important point being that Mi = lOTeV is within the acceptable range. 

We now consider the raw asymmetry one order more than the previous case i. e. 7f™ L = 
a x 10~ 9 . From fig. 14.31 we see that for rh\ = 10~ 3 eV there is only an upper bound 
Mi = lO 10 ' 5 ^^, such that the final asymmetry lies in the allowed region for all smaller 
values of Mi. Thus for the case of raw asymmetry an order of magnitude smaller, the upper 
bound on Mi decrease by two orders of magnitude (e.g. compare previous paragraph). 
However, the choice of smaller values of rh\ leads to a small window for values of Mi for 
which we end up with the final required asymmetry. In particular for rh\ = 10~ A eV the 
allowed range for Mi is (10 12 — 10 13 )Ge^, while for fh\ = 10~ 5 eV the allowed range 
shifts to (10 14 — 10 15 )Gel^. The window effect can be understood as follows. Increasing 
the value of M 1 tends to lift the suppression imposed by the m 2 dependence of the wash 
out effects, thus improving efficiency of the latter. However, further increase in Mi makes 
the effects too efficient, erasing the raw asymmetry to insignificant levels. 
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Figure 4.3: The allowed values of Mi against the final required asymmetry is shown for 
7] raw = 2.0 x 1(T 9 

The windowing effect emerges clearly as we consider the cases of large raw asymmetries. 
This is shown in fig. 14.41 It is seen that as the raw asymmetry increases the allowed regions 
become progressively narrower and lie in the range (10 10 — lO 15 )^^. Thus a given raw 
lepton asymmetry determines a corresponding small range of the heavy Majorana neutrino 
masses for which we can obtain the final asymmetry of the required order a x 10~ 10 . Again 
smaller is the effective neutrino mass rhi larger is the mean value of the allowed mass of 
the heavy Majorana neutrino and this is a consequence of normal see-saw. 

Finally, in the following, we give an example for non-thermal creation of L-asymmetry 
in the context of left-right symmetric model. 

4.4 Lepton asymmetry in left-right symmetric model 

We discuss qualitatively the possibility of lepton asymmetry during the left-right symme- 
try breaking phase transition . In the following we recapitulate the important aspects 
of left-right symmetric model for the present purpose and the possible non-thermal mech- 
anism of producing raw lepton asymmetry. This asymmetry which gets converted to 
baryon asymmetry, can be naturally small if the quartic couplings of the theory are small. 
Smallness of zero-temperature CP phase is not essential for this mechanism to provide 
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Figure 4.4: The allowed region of Mi is shown for different values of m\ for large values 

of r] raw 

small raw L asymmetry. 

4.4.1 Left-Right symmetric model and transient domain walls 

The important features of the left-right symmetric model based on the gauge group 
SU{2)l x SU{2)r x U{1)b-l are elucidated in chapter 3. The Higgs potential of the 
theory naturally entails a vacuum structure wherein at the first stage of symmetry break- 
ing, either one of A^ or A# acquires a vacuum expectation value the left-right symmetry, 
SU(2)l <-> SU(2)r, breaks. It is required that Ar acquires a VEV first, resulting in 
SU{2) R <g> U(l) B -L -> U{l) Y . Finally Q = T\ + T| + \{B - L), survives after the 
bidoublet $ and the acquire VEVs. 

If the left-right symmetry were exact, the first stage of breaking gives rise to stable 
domain walls jHOl EH E2] interpolating between the L and R-like regions. By L-like we 
mean regions favored by the observed phenomenology, while in the R-like regions the 
vacuum expectation value of A# is zero. Unless some non-trivial mechanism prevents 
this domain structure, the existence of R-like domains would disagree with low energy 
phenomenology. Furthermore, the domain walls would quickly come to dominate the 
energy density of the Universe. Thus in this theory a departure from exact symmetry in 
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such a way as to eliminate the R-like regions is essential. 

The domain walls formed can be transient if there exists a slight deviation from exact 
discrete symmetry. As a result the thermal perturbative corrections to the Higgs field free 
energy will not be symmetric and the domain walls will be unstable. This is possible if the 
low energy (~ 10 4 GeV-10 9 GeV) left-right symmetric theory is descended from a Grand 
Unified Theory (GUT) and the effect is small, suppressed by the GUT mass scale. In the 
process of cooling the Universe would first pass through the phase transition where this 
approximate symmetry breaks. The slight difference in free energy between the two types 
of regions causes a pressure difference across the walls, converting all the R-like regions 
to L-like regions. Details of this dynamics can be found in ref. |76j . 

4.4.2 Leptogenesis mechanism 

In order to produce adequate lepton asymmetry the following criteria of Skharov's have to 
be satisfied [T]. C and CP- violation and finally all the thermal processes have to be out of 
equilibrium. The first two properties can be realized if there is a CP- violating condensate 
exists in the domain wall. Finally the out-of-equilibrium condition can be realized from 
the directional motion of the domain walls to R-like regions and thus making all the 
domains L-like. 

We now consider the interaction of neutrinos from the L-R wall, which is moving 
towards the energetically disfavored phase, the R-like region. The left-handed neutrino, 
z/^s, are massive in this domain as they couple with their CP-conjugate states, where as 
they are massless in the phase behind the wall as (Ax,) = 0. This can be seen from the 
Yukawa coupling 

C yuk = fA L u c L u L + h.c. (4.14) 

and since Al has a wall like profile, the mass of is zero behind the wall and is 0{vr) 
in front of it. 

To get a net lepton-asymmetry one needs the asymmetry between the reflection and 
transmission coefficients from the wall between and its CP-conjugate state (^£). If it 
favors the transmission of to the L-like region then the excess of antineutrinos (z/£) 
reflected in front of the wall will equilibrate with the v L due to helicity flipping scatterings. 

49 



However the transmitted excess of Vl remain as it is, since it does not couple with v c L 
behind the wall. 

At least two of the Higgs expectation values in L-R model are generically complex, 
thus providing natural CP-violation [H3] permitting all parameters in the Higgs potential 
to be real. We choose the following profiles for the Higgses 

(Al)= I I A /e ). | ° °| (1.10) 





v R / 

In the classical approximation where the wall width is assumed to be large the CP violat- 
ing phase inside the domain wall becomes position dependent. Under these circumstances 
a formalism exists [H3J IE3J EH] , wherein the chemical potential fi L created for the Lepton 
number can be computed as a solution of the diffusion equation 

-D v \i' L - v w fi' L + d(x) T hf /i L = S(x). (4.16) 

Here D u is the neutrino diffusion coefficient, v w is the velocity of the wall, taken to be 
moving in the +x direction, Thf is the rate of helicity flipping interactions taking place in 
front of the wall (hence the step function 0(x)), and S is the source term which contains 
derivatives of the position dependent complex Dirac mass. 

After integration of the above equation and using inputs from the numerical solutions 
we find the raw Lepton asymmetry [ZE] 

1 M 4 

raw ^ 0.01 V w — -=r±- (4.17) 

/L w g, T 5 A W 1 ; 

where 7]™ w is the ratio of to the entropy density, s = (2ir 2 / '45) g*T 3 . In the right 
hand side Mi stands for the Majorana neutrino mass, A w is the wall width and g* the 
effective thermodynamic degrees of freedom at the epoch with temperature T. However, 
the high temperature sphalerons are efficiently converting the L asymmetry into B — L 
asymmetry. The standard chemical equilibrium calculation 03] gives rf™ L = t^t]™ w . 
Using Mi = fA T , with A T is the temperature dependent VEV acquired by the A R in the 



phase of interest, and A m ! = ^X e ffA T in equation (j4.17J) we get 



— ~10-V(^) 5 Z 4 ^- (4.18) 
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Here we have used g* = 110. Depending on the Majorana Yukawa coupling the raw 
asymmetry can take a range a values of 0(1CT 4 — 1CT 10 ). 

4.5 Results and discussions 

We have assumed a non-thermal production of raw lepton asymmetry during the B — L 
breaking phase transition in a generic theory incorporating B — L as a gauge symme- 
try. If this asymmetry passes without much dilution to be the currently observed baryon 
asymmetry consistent with WMAP and Big Bang nucleosynthesis, then the effective neu- 
trino mass parameter rhi must be less than 10~ 3 eV. Solution of the relevant Boltzmann 
equations shows that for rhi = 10~ 4 eV the mass of lightest right handed neutrino Ni has 
to be smaller than 10 12 GeV and can be as low as 10 TeV. In a more restrictive scenario 
where the neutrino Dirac mass matrix is identified with that of the charged leptons it is 
necessary that Mi > 10 7 GeV in order to satisfy rhi < !0~ 3 eV. Therefore in the more 
restricted scenario all values Mi, 10 8 GeV < Mi < 10 12 GeV can successfully create the 
required asymmetry. If the Dirac mass scale of neutrinos is less restricted, much lower 
values of M x are allowed. In particular, a right handed neutrino as low as 10 TeV is 
admissible. 

If the raw asymmetry is large, the numerical solutions show a small window for M\ to 
get the final asymmetry of the required order. The allowed range gets smaller as the raw 
asymmetry gets larger. This is true for all allowed values of the neutrino mass parameter 
rhi. 

In summary, if the B — L symmetry is gauged, we start with a clean slate foiB — L 
number and an asymmetry in it can be generated by a non-perturbative mechanism at 
the scale where it breaks. The presence of heavy right handed neutrinos still permits 
sufficient asymmetry to be left over in the form of baryons for a large range of values of 
the B — L breaking scales. While other mechanisms of leptogenesis become unnatural 
below 10 7 GeV (see for example chapter 1 and chapter 2) this mechanism even tolerates 
TeV scale. A specific mechanism of this kind is possible in the context of Left-Right 
symmetric model, presumably embedded in the larger unifying group 50(10). Here we 



51 



conjecture that upon incorporation of supersymmetry, the qualitative picture will remain 
unaltered and the present gravitino bound of 10 9 GeV for reheating temperature after 
inflation can be easily accommodated. 
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Part-II 



Topological Defects 
Cosmic Strings 
and 

Baryogenesis via Leptogenesis 
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Chapter 5 

Soliton fermion systems in Quantum 
Field Theory 

5.1 Introduction 

"Solitary waves" and "solitons" emerge as the non perturbative solutions of non-linear 
wave equations in classical field theory. These are non-dispersive localized packets of 
energy moving uniformly in space. Although they have many common features, the former 
class of solutions do not keep their shape intact during a collision among themselves, 
whereas the latter do in the asymptotic time domain. Thus the latter class of solutions is 
a special subset of the former class of solutions, but the converse is not necessarily true. 
Having differentiated solitons from solitary waves, we should mention (without embarking 
into details) that the distinction between them in field theory is completely blurred. So 
here onwards, what merely called as "solitary wave" , we call it "soliton" as it is more crisp 
and appealing. A prototype example in field theory is the well known massive Thirring 
model where the solution is exactly solitonic. 

It was believed that elementary particles in nature can be thought of as localized 
packets of energy, the prime among such theories is the famous Skyrme model. But 
till date no satisfactory formalism has been developed to view elementary particles as 
solitons as the present effective theory takes particles as point like, nor is there any 
experimental evidence at latest attainable energies to prove that particles have indeed 
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extended features. We bypass this issue here but concentrate on defects occurring in gauge 
theories. Of course field theories describing elementary particles are quantum theories, 
whereas solitons are solutions of classical field theories. Although it is a classical solution, 
the quantization of bosonic fluctuation in the background of this non trivial vacuum, the 
so called "semiclassical treatment" , leads to a finite correction [SH] • 

A topological solution can be obtained as a solution of a differential equation with 
a difficulty that the boundary values are known only at infinity. In the context of field 
theories a prototype example is the 'kink', which emerge as the solution of 4 -theory in 
1+1 dimensions. For demonstration purpose, in section 5.2 we consider a toy model in 
1 + 1 dimensions. 

A more curious phenomenon happens in Quantum Field Theory (QFT) when a fermion 
interacts with a kink. This was first pointed out by Jackiw and Rebbi in 1976. It 
was shown that the resulting fermion number becomes fractionalized. Although it is 
surprising, the consequences of this phenomena have already been studied extensively in 
literature in the context of condensed matter systems, polyacetylene being the standard 
example [EH EH IH3 EH] , and in high energy physics, cosmic string being the example [SUl 



5.2 Solitary waves in 1+1 dimensions 

Soliton or solitary wave is a static, localized and finite energy solution of non-linear wave 
equation in classical field theory. For the present purpose we consider only bosonic fields 
<f)(x,t) in one space and one time dimensions. The dynamics of the fields is governed by 
the Lagrangian 



EH- 



C = ~d fl <f>d»<f>-V( ( f)), 




where V{4>) is the required potential. In A0 4 -theory it is g: 



iven by 





The Hamiltonian corresponding to the Lagrangian ()5.1|) is 



n= l -d,ct>d^ + v{ ( t>). 



(5.3) 
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In order to get solitonic solutions we demand that is independent of time. Thus the 
equation of motion will be 



d 2 _ dV{(f>) 
dx 2 dip 

Again the fmiteness of the energy requires that the energy functional 



(5.4) 



E[<p} = J^°°dx U^j + V{<j>) <oo. (5.5) 

If the absolute minima of V(4>), which are also its zeros, occur at n points i.e 

V(4>) = for <p = g\i = 1,2 ---n (5.6) 

then the energy functional E[<f>] can be minimized when the field </>(x,t) is constant in 
space-time and takes any one of these values. That is E[(/>] = if and only if <p(x, t) = g l . 
Further to make E[<p} finite at spatial infinity (x — > ±oo), it requires that V{<f>) = at 
x — > ±oo, which is the absolute minimum of V(<f>). In addition to that ^ = 0. This 
implies that 0=constant. 

We now solve Eqn. (j5.4|) . One of the easiest way to solve this equation is that, if we 
replace 

4> «-> y 

X <-> t 
V{cj>) ~ -V{cj>) 

in equation (|5.4j) then it will resemble with newtons second law of motion of a particle 
of unit mass moving in a negative potential. The solution <p{x) represents this particle's 
motion of this analogue particle. So the total energy of this motion is conserved as x, the 
time, varies and is given by 

i '•■-■> 



W =2[t) (57) 
Note that this energy is different from the energy of the soliton. Now using the above 
boundary conditions (used for soliton), we have W = 0. This gives rise to 

/ dx = ± —f==- 5.8 

Jo Jm y f W{^) 
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Solving this equation for the particular potential defined in Eqn. (|5.2|) we get 

4>±( x ) — ±— 1= tanh(//a;). (5.9) 
v A 

The solution corresponding to <fi+(x) is called 'kink' and that of <p-{x) is called 'antikink'. 
Note that the coupling constant A appear in the denominator of the above Eqn. (j5.9|) . If 
A — > then the corresponding solution — > oo. Thus the phenomenon is completely 
non-perturbative one. The energy corresponding to these solutions is given by 



E= edx, (5.10) 

J — oo 

where 

e = — sec/i 4 (yux). (5-H) 
A 

Sometimes this energy is called classical 'kink mass' M c \ = E = At this point it is 
worth to point out the symmetric properties between kink and antikink. It is clear from 
equation ()5.9|) that 

</>+(x) = -<f>-(x) 

= <p-{-x). (5.12) 

Thus a kink starts its motion from a trivial vacuum <p = — ^ at i ^ - ooto reach 
other trivial vacuum <ft = at x — > +oo through zero. Although these are extended 
objects still its energy is concentrated within narrow region of space of 0(1/ fi). So these 
solutions is expected to behave as particles in high energy physics. Fortunately these 
solutions admit a Lorentz boost with a definite velocity (v). For example, the energy of 
the soliton transforms as 

E - E' = -j^L=. (5.13) 

V 1 — IT 

It is worth to point out that kink is merely a solitary wave and not a soliton. That 
means they do not survive collisions. Since in our current discussion (i.e. in context of 
4 -theory) there is only a pair of kink and antikink solution exist, so it is difficult to check 
their shape retaining property in a collision. Of course, a numerical calculation is needed 
to check it. But it is beyond the scope of this manuscript. 
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5.3 Topological charge 



Often we need to distinguish the topological solutions of a system of equations. So we 
define a 'topological index', which is conserved in time. It is very similar to 'quantum 
number' in QFT, but its origin is completely different. Now we will illustrate it in context 
of A0 4 -theory. The potential has two degenerate minima, at <fi = Consequently 
solutions of this system, whether static or time dependent, fall into four topological sec- 
tors. These are characterized by the pair of indices (— -4=, -4=), (-7=, — -7=), (— -7=, — ^), 
( v^x' res P ec ti ver y- Where the quantities inside the parenthesis represents the value of 
4>±{x) at their spatial infinity. Now we choose a particular sector (— — -^=), which is 
nothing but the value of kink at x — > — 00 and antikink at x — > 00. Even though we may 
not be able to calculate easily what happens after they collide, but we can be sure that 
the resulting field configuration will lie in the (— -^=, — sector. At this juncture it is 
worth to designate such sectors by a definite quantity called 'topological charge', defined 
by 

\/A 

Q = J_ = +oo) - 4>(x = -00)] , (5.14) 



associated with a conserved current 



k a = y^ e ^d^ (5.15) 



where a, (3= 0, 1. From equation (|5.15J) we get the topological charge 



+00 



Q = / k dx. (5.16) 
J —00 

Thus we need both 0(oo) and </>(— 00) to identify a topological sector. If Q = 0, then the 
object is a non topological object. Thus topological charges are nothing but the boundary 
conditions of the problem in contrast to usual Noether charge in QFT which comes from 
the continuous symmetry associated with the theory. 

5.4 Fractional fermion charge in QFT 

The phenomenon of fractional charge was first discovered by Jackiw and Rebbi [17] in 
their pioneering work in 1976. We begin by reviewing their work. Note that the charges 
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we are talking about are not the wellknown and widely accepted fractional charge of 
quarks. But it is the charge which appear as the eigen value of number operator in the 
QFT. For our illustration purpose, we consider the 1+1-dimensional theory. 

As an extension of equation (|5.1jl to include fermion field if), the additional terms are 
given by the Lagrangian 

£fermion = 4> (*Y '% ~ 9<j>fat)) i> • (5-17) 

For the potential defined in equation ([5.2)1 . the trivial vacuua are <f> = It is well 

known that the presence of two degenerate classical solutions indicate spontaneous break- 
ing of the — > —(f) symmetry of the Lagrangian ()5.1jl . Around each of these classical 
vacuua a whole tower of Fock states can be built and hence to be called vacuum sectors. 
In addition to that the system has also two other static solutions. Those are the solitonic 
solutions, the so called kink and its reflected partner the anti-kink, given respectively by 

S = tanh(yux). (5.18) 

As per the general theory of semi classical quantization of quantum fields, one can build 
two other separate towers of states, one around each of these solitonic solutions. In short, 
we have four sectors of states for this interacting system: two are the vacuum sectors built 
around (f> = ±-7= and others two are the soliton sectors ±(/x/\/A) tanh(/ix). 

5.4.1 The vacuum sector 

In the vacuum sector, <f> = the Lagrangian for the fermion field if) is given by 

C, = ^(i 1 v d v -m F )if), (5.19) 

where m F = g(fx/ a/A) is the mass of the fermion and v stands for 0,1. Since we have only 
one space dimension, we use the representation of the Dirac matrices as 7 1 = (3a = io~ 3 
and 7 = (3 = cr 1 . Let us denote by Uk{x) and Vk(x) the positive and negative energy 
spinorial solutions of the Dirac equations 

(-iad x + (3m F )u k {x) = E k u k (x) (5.20) 
{-iad x + (3m F )v k (x) = -E k v k (x), (5.21) 
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where E k = +\Jk 2 + m? F and spinor indices have been suppressed. The Dirac field can 
be expanded in terms of these solutions and the destruction operators b k and d k , obeying 
the usual anticommutation rules, as 

i/>{x,t) = Yy> k u k e~ iEkt + d\v k e lEkt }. (5.22) 
k 

The vacuum state in the <fi = (+/i/a/A) sector is given by the conditions 

6*10) = 4|0) = (5.23) 

with all the bosonic oscillators being in the ground state. 

Note that the third Pauli matrix a 3 acts as the charge conjugation matrix. It anticom- 
mutes with the Dirac Hamiltonian in Eq. (j5.21jl and generates for every positive energy 
solution Uk(x) of energy Ef. the corresponding negative energy solution Vk{x) of energy 
-E k , i.e. 

a 3 u k (x) = v k {x). (5.24) 

Hence all modes of the expansion ()5.22j) come in pairs with equal positive and negative 
energies. There are no zero energy solutions for tp in the vacuum sector that allowed by 
the Dirac equations (|5.20j) and f)5.21D . 

Finally we consider the number density operator 

p{x,t) = ^[^{x,t),tfj{x,t)}. (5.25) 

The form of p is designed so as to be odd under charge conjugation. Inserting the mode 
expansion ()5.22j) and using the orthonormality conditions of the Dirac spinors, the total 
charge becomes 

Q = fd M x,t) 

= iT.iV'lM + I*.4l) 

k 

= E (\f>l b * - X /2] - [44 - 1/2]) 

k 

= (5-26) 
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Notice that the half-integers cancel term by term because of the existence of paired positive 
and negative energy modes. Hence the familiar result in the vacuum sector that the charge 
operator has only integer eigenvalues. 



5.4.2 The soliton Sector 

In this section, we repeat the same procedure as we did in the trivial vacuum sector. We 
replace the by its solitonic value cf> s from Eq. (j5.18|) . The corresponding Dirac equation 
becomes 

(— iad x + (3m F tanh(/ix)) ip(x) = Eip(x). (5-27) 

This equation has a set of positive energy solutions ipk( x ) associated with energy E k . 
However, the charge conjugation matrix a 3 again anticommutes with the Dirac Hamil- 
tonian, C~ l HC = —H. Therefore, for every positive energy solution ipk( x ) there is a 
negative energy solution ipk{x) with energy — Ek- 

For the two components ipi,2 of the Dirac spinor, equation (|5.27|) yields the coupled 
equations 

{—d x + m F tanh(fix)) if} 2 = Eip 1 (5.28) 
(d x + m F tanh(/ix)) ipi = Eip 2 - (5.29) 

Solving these coupled equations we get an unpaired zero-energy solution 

_ I Nex P {~ m F Jo dx' tanh(/ix')) J ^ 

\ /' 

In contrast to the normal vacuum, in the soliton sector we have a zero energy solution 
([5.30)1 . This is because of the soliton function (fi/y/X) tanh(//x) which forms the back- 
ground potential for the Dirac spinor tends to opposing limits ±(/i/v / A) as x — > ±oo. 
In infinite spatial volume this solution has no partner. It is self charge conjugate i.e. 
^i^o — "0Q- The mode expansion of the Dirac field operator now becomes 

V(x,t) = Yl [hVk(x)e- tEkt + 4%(x)e ii ^] + a^{x) (5.31) 

where 'a' is the destruction operator for the zero mode ip . 
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Unlike the vacuum sector built around (/x/a/A) which had a unique ground state, in 
the solitonic sector there exists two degenerate ground states because of the existence of 
a zero energy solution of ip. They are |so/)_, being the unfilled state, and \sol) + , being 
the filled state, obeying 

afc|so/)_ = bk\sol)- = dk\sol)^ = (5.32) 

and 

\sol)+ = 4|so/)_ 
a\sol)+ = \sol)-. (5.33) 

These are the two basic quantum soliton states of this system. They are energetically 
degenerate, but are distinguishable by their charge 

Q = \J dx[^(x,t),i;(x,t)] 

= ^([ & I'^ + [4,4]) + (l/2)[a t ,a] 

k 

= E - V2) - (44 - 1/2)) + (4a - 1/2) 

k 

= E fa* ~ 44) + (4a - 1/2) (5.34) 

k 

Notice that the piece (-1/2) coming from the zero mode commutator remains uncanceled 
because it does not have a charge conjugate partner. Thus it is obvious that the total 
charge (the number operator) Q has half integral eigen values. The two degenerate soliton 
states have eigen values ±1/2 respectively for the total number operator Q: 

Q\sol)_ = -(l/2)|soZ>_ 

Q\sol)+ = (l/2)\sol)+ (5.35) 

Thus the solitonic states are superselected from the normal vacuum due to the unavail- 
ability of fractional states in the translational invariant vacuum. Note that the fraction- 
alization of the fermion number in the background of a nontrivial vacuum is one of the 
manifestation of Dirac negative energy sea. 
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5.5 Implication of fermion zero modes and superse- 
lection rules 



By the name 'superselection' we mean that there are certain restrictions on nature and 
scope of possible measurements. It is shown that [92J there exist superselection rules for 
spinor fields which can be proved in context of the assignment of parity to the quantum 
mechanical states of a system. 

The usual assumption in quantum mechanics is that it is possible to carry out a 
complete set of measurements; the result of which determines the state of any vector 
completely except for the usual phase factor. If there exist sectors A, B, C, etc. in 
the Hilbert space such that vectors in each subspace can be independently rescaled by 
phases, such sectors must be understood to be completely independent, and superposition 
of states from such unrelated sectors should be forbidden. Specifically, let state vectors 
^/a, ^b,^c etc. belong to such independent sectors, so that it is assumed that no physical 
measurements distinguish between the state 

*a + *b + *c + --- (5-36) 

and 

e ia ^ A + e^^B + e^^c + ■■■ (5.37) 

where a, (5, 7, . . . are arbitrary phases. Then the expectation value of any operator 
possessing matrix elements connecting subspaces A and B, or A and C, etc. must be 
completely undefined. Hence such an operator will not correspond to a measurable quan- 
tity. 

It is customary to say that a selection rule operates between subspace of the total 
Hilbert space if the state vectors of each subspace remain orthogonal to all state vectors 
of other subspaces as long as the system is isolated. There is, for instance, a selection 
rule which prevents any state of an isolated system from changing its total linear momen- 
tum. Similarly, the state vectors of the subspace containing all states with total angular 
momentum J will remain, in a closed systems, orthogonal to all states with any other 
total angular momentum. So we shall say that a superselection rule operates between 
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subspaces if there does not exist any selection rule between them and if, in addition to 
this, there are no measurable quantities with finite matrix elements between their state 
vectors. For example if we superpose states of integral spin (bosonic) with half integral 
spin (fermionic) then the probability of finding of each states are equal and 1/2 each. 
So it is impossible to compare directly a spin half particle with a spin one particle. It 
is conjectured that [92J superselection rule not only work between the states of different 
intrinsic parity but also it works between the states of different total charge. For instance 
if we construct a state, by superposing the states |l/2) and |1), as 



then it is expected that under P 4 operation it will come back to it's original state up to a 
phase. To verify this we can assign a parity ±1 to a Dirac neutrino and ±i to a Majorana 
neutrino in the trivial vacuum sector. But in the vortex sector the corresponding state is 
doubly degenerate and thus the parity assignment is half of their trivial vacuum values. 
Now the P 4 operation on (|5.38|) for a Majorana neutrino will give 



This state is another state orthogonal to ()5.38|) and hence implies the impossibility of 
superposition of such states. We discuss more about it in context of the decay of cosmic 
strings in chapter 6. 



^ = -[|l/2) + |l>] 




PV. = ^[-|i/2> + |i>]. 



(5.39) 
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Chapter 6 

Metastable topological defects and 
fermion zero modes 

6.1 Introduction 

In the process of cooling down from the highest possible energy scale till the present epoch, 
early Universe has been passed through several phase transitions. A phase transition 
occurs through the spontaneous breaking of a larger gauge symmetry group G to any of 
its subgroup H. If the corresponding vacuum manifold M. = G/H is non trivial then the 
topological defects are formed. Depending on the geometry of the manifold the possible 
defects are domain walls, cosmic strings, monopoles and textures. Domain walls form 
if M. has disconnected components, strings can form if M. has unshrinkable loops and 
monopoles form when Ai contains unshrinkable surfaces. The relevant properties of the 
manifold M. are most conveniently studied using homotopy theory; the n th homotopy 
group ir n (AA) classifies qualitatively distinct mappings from the n-dimensional sphere S n 
into the manifold M. If n (M) ^ I then the formed defects are called domain walls. 
On the other hand, if it\(M) ^ I then cosmic strings are formed and if ir 2 (M) ^ I then 
monopoles are formed. These are extended objects, cosmic strings being 1+1-dimensional, 
domain walls the 2+1-dimensional and finally monopoles the 3+1-dimensional but not 
point like. 

The topological defects emerge as the solitonic solutions of non linear wave equations in 
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gauge theories. Therefore, in the asymptotic time domain these objects are very stable and 
keep their original configurations in tact. However, they can be unwound by overcoming 
finite energy barriers. The dynamics of these objects may have several consequences that 
may link to the physics of early Universe with the recent cosmology. In particular, we 
discuss baryogenesis via leptogenesis from the decay products of cosmic strings in chapter 
7. 

6.2 Metastable topological defects 

If there are several stages of symmetry breaking, the topological stability of these ex- 
tended objects depend on the structure of the vacuum manifold at every stage. Under 
certain conditions defects stable at the first stage of symmetry breaking are rendered 
unstable at lower temperatures. Likewise objects that seem to enjoy topological stability 
in the low energy effective theory are actually unstable in the complete description of 
the theory. The instability however tends to energetically exorbitant thus rendering the 
object metastable 

There are several reasons that may stabilize objects that are unstable on topological 
grounds. For instance superconducting strings lead to small loops called vortons, stabi- 
lized by the electric current flowing through them. Here we consider the occurrence of 
fermionic zero-energy modes. As a rule of thumb, if the number of trapped zero-energy 
modes is n, the ground state of the defect carries fermion number n/2. Thus if n is odd 
one obtains curious occurrence of fractionalization of fermion number, proven to be inte- 
ger in translation invariant field theory. In general there are several internal charges the 
fermions may carry and the fractionalization rule affects all of them. It has been observed 
by several authors [HU |H3 EH] that the unavailability of final states of matching charge 
should forbid the decay of such objects in isolation. 

Here we consider the sequential breakdown of the gauge group G — > Hi — > H 2 - If 
the vacuum manifold M. = G/H\ is non-trivial then topological defects are formed. 
However, these defects in their corresponding low energy theories loss their stability after 
the breaking of second phase of symmetry if the manifold G/H 2 is simply connected and 



66 



thus making the defect metastable. As an example we study the metastable cosmic strings 
in sec. 6.2.1. 

6.2.1 Metastable cosmic strings and fermion zero modes 

Here we construct two examples in which the cosmic strings of a low energy theory are 
metastable due to the embedding of the low energy symmetry group in a larger symmetry 
group at higher energy. Examples of this kind were considered in [§3]. Borrowing the 
strategies for bosonic sector from there, we include appropriate fermionic content to ensure 
the zero-modes. 

A. Dirac fermions 

Consider first a model with two stages of symmetry breaking similar to j^H], but with 
local £77(2) gauge invariance. The two scalars £ and a are respectively real triplet and 
complex doublet. The Lagrangian is taken to be 

£ = -~F^ a F£ v + ~Dj:- D^ + D^D^a 
- A 1 (S-S-r /l 2 ) 2 -A 2 (ata-r ?2 2 ) 2 

+ Ai 2 ^iE • (Tto (6.1) 

where isovector notation is used and r are the Pauli matrices. It is assumed that r\\ ^> i]2 
and that the coupling Ai 2 > satisfies \12V1 ^ ^2^2 

The vacuum expectation value (VEV) E = (0 t]i) T breaks the SU(2) to the £7(1) 
generated by exp(2T 3 a/2). The effective theory of the a can be rewritten as the theory 
of two complex scalar a u and ad for the up and the down components respectively. The 
potential of the effective theory favors the minimum 

cr u = &d = (6.2) 

In the effective theory a u enjoys a £7(1) invariance a u — *> e ia a u which is broken by the 
above VEV to Z = {e 2nnl }, n = 1,2, . . .. This makes possible vortex solutions with an 
ansatz in the lowest winding number sector 

a u (r,4>)=r ]2 f(r)e- i t> (6.3) 
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where r, <p are planar coordinates with vortex aligned along the z axis. The vortex 
configuration is a local minimum, however it can decay by spontaneous formation of 
a monopole-antimonopole pair [OH]- These monopoles are permitted by the first break- 
ing SO (3) — > 5*0(2) in the £ sector. Paraphrasing the discussion of (HS], we have 
SU(2) — >U(1) — > Z. The vortices are stable in the low energy theory because 7r 1 ([/(l)/Z) 
is nontrivial. But in the SU(2), the Z lifts to { e 4n7rlr3 / 2 } = / making it possible to unwind 
the vortex by crossing an energy barrier. 

Consider now the introduction of a doublet of fermion species = (N L , E L ) assumed 
to be left handed and a singlet right handed species Nr. The Yukawa coupling of these 
to the a u is given by IiNr/t^l, which in the vortex sector reads 

Cc-i> ~ /^ 2 /(r)(e-^AQV L + h.c.) (6.4) 

The lowest energy bound states resulting from this coupling are characterized by a topo- 
logical index, j3T] X = ni — hr where ni and hr are the zero modes of the left handed 
and the right handed fermions respectively. This index can be computed using the for- 
mula inn ezi 

I=-L(lndetM)|^ (6.5) 
2m r 

where M is the position dependent effective mass matrix for the fermions. In the present 
case this gives rise to a single zero-energy mode for the fermions of species N. According 
to well known reasoning jUj to be recapitulated below, this requires the assignment of 
either of the values ±1/2 to the fermion number of this configuration. 

B. Majorana fermions 

The example above can be extended to the case where the Nr is a Majorana fermion. 



Being a singlet N R admits a mass term M m N^Nr, Mm signifying Majorana mass. This 
could also be a spontaneously generated mass due to the presence of a neutral scalar x 



with coupling terms HmxNrNr + h.c If this x acquires a VEV at energies higher than 
the S, the ./V^ particles possess a Majorana mass and fermion number is not a conserved 
observable. 
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Finally we present the case where Majorana mass is spontaneously generated at the 
same scale at which the vortex forms. Consider a theory with local SU(3) symmetry 
broken to U(l) by two scalars, $ an octet acquiring a VEV 771 A3 (A3 here being the third 
Gell-Mann matrix) and 0, a 3, acquiring the VEV (<p k ) = r)2<5 fc2 , with 772 <C f]i- Thus 

SU(3)-^U(1) 3 ® U(1) 8 ^U(1) + (6.6) 

Here U(l) 3 and Z7(l) 8 are generated by A 3 and A 8 respectively, and U(l) + is generated 
by (v^As + Aa)/2 and likewise ZV( 1)_ to be used below. It can be checked that this 
pattern of VEVs can be generically obtained from the quartic scalar potential of the 
above Higgses. The effective theory at the second breaking U(l)- — ► Z gives rise to 
cosmic strings. However the Z lifts to identity in the SU(3) so that the string can break 
with the formation of monopole-antimonopole pair. 

Now add a multiplet of left-handed fermions belonging to 15. Its mass terms arise 
from the following coupling to the 3 

^Majorana = VfV for^O (6-7) 

The indices symmetric under exchange have been indicated by curly brackets. No mass 
terms result from the 8 because it cannot provide a singlet from tensor product with 
15 <S> 15 PHI- After substituting the <fi VEV a systematic enumeration shows that all but 
the two components ip{ 22 ^ and ip^ 22 ^ acquire Majorana masses at the second stage of the 
breaking. Specifically we find the Majorana mass matrix to be indeed rank 13. Thus, 
using either of the results [SU] or i.e., eq. ([6.5|) we can see that there will be 13 zero 
modes present in the lowest winding sector of the cosmic string. Thus the induced fermion 
number differs from that of the vacuum by half-integer as required. 

C. Final state zero-modes 

The stability argument being advanced is in jeopardy if the final state after rupture of the 
topological object also possesses half-integral fermionic charge. To see that this is not the 
case it is necessary to study the zero-modes on the two semi-infinite strings shown in fig. 
16.11 Generically we expect each of the halves to support the same number of zero-modes, 
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making the total fermion number of the putative final state integer valued, as required 
for the validity of our argument. 



z 
A 



Figure 6.1: Schematic configuration of isospin vectors after the rupture of a string. Inter- 
nal orientations are mapped to external space. They are shown just outside the core of 
the two resulting pieces and on the mid-plane symmetrically separating the two. 

Consider the ansatz for the lower piece (/) with origin at the corresponding monopole 
and coordinates (r h 9i,<f)). For the domain z < let the ansatz for the field a be 



ur 



° //(/V) exp{>7'-0} ( ° | /,(,-,) 



(6i 



so that {a) has the behavior 



fiin) for e l « o 



,12, 



f]2e 



-icf> 







fi(ri) for 9i w vr 



which agrees with the ansatz (J6.3)) for the cosmic string at the South pole. Likewise for 
the domain z > 0, i.e. the upper piece (u), we choose 



U 



(6.9) 
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resulting in the behavior 




f u (r u ) for 9 U « 




for 0,, 



7T 



thus matching correctly with the cosmic string at the North pole. The ansatz for the 
heavier scalar £ needs to be appropriately set up, C/°°S ■ tU°°^ in both / and u domains. 
This scalar however does not contribute to fermion mass matrix. 



The two maps match at the mid-plane where Q[ = n—9 u and a u ~ e at (9; = 9. 



7T 



/2 



so that we have ensured that the combined map is within the same homotopy class as 
the string we began with. Finally, as the two pieces move far away, each can be seen to 
have the same number of zero-modes. To see this we can choose [HBJ EH] fermion ansatz 
for the zero-modes compatible with the scalar field ansatz, in each of the patches I and u. 
In (isospin)®(two component spinor) notation for ip£ and for the two component fermion 
N R , 



u c 




n. 



R 




(6.10) 



where the labels I, u have been dropped. To analyse the asymptotic radial dependence 
choose 7 r = a 2 the Pauli matrix. In each patch one finds the pair y?i(r), X2{f) ~ e~ hri2r to 
constitute the zero-mode while for the other pair, ^{f), Xi( r ) ~ e +hr,2r which are therefore 
not normalizable. In any case, since each of the pieces acquires the same number of zero- 
modes, the total fermion number of the putative final state has been proved to be integer 
as required. 



6.3 Assignment of fermion number 

We now recapitulate the reasoning behind the assignment of fractional fermion number. 
We focus on the Majorana fermion case, which is more nettlesome, while the treatment 
of the Dirac case is standard [571 HOOj . In the prime example jSOl in 3 + 1 dimensions of 
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a single left-handed fermion species \I/ l coupled to an abelian Higgs model according to 

C i) = iW~ Ll »D^ L - l -{ha¥lm L + h.c.) (6.11) 

the following result has been obtained. For a vortex oriented along the z-axis, and in the 
winding number sector n, the fermion zero-modes are of the form 

^ L0 (x) = I 1 1 [I7(r)e fl * + lA*( r ) e i(»-i-0^1 gi ( z + t ) ( 6 .12) 

w 

In the presence of the vortex, r 3 (here representing Lorentz transformations on spinors) 
acts as the matrix which exchanges solutions of positive frequency with those of negative 
frequency. It is therefore identified as the "particle conjugation" operator. In the above 
ansatz, the \&£ in the zero-frequency sector are charge self-conjugates, t 3 ^ L0 = ^lo, an d 
have an associated left moving zero mode along the vortex. The functions satisfying 
t3 ^lo — ~^lo are not normalizable. The situation is reversed when the winding sense 
of the scalar field is reversed, ie, for <j u ~ e~ in ^. In the winding number sector n, regular 
normalizable solutions jSU] exist for for < I < n — 1. The lowest energy sector of 
the vortex is now 2 n -fold degenerate, and each zero-energy mode needs to be interpreted 
as contributing a value ±1/2 to the total fermion number of the individual states |q"U] . 
This conclusion is difficult to circumvent if the particle spectrum is to reflect the charge 
conjugation symmetry of the theory |lUlj . The lowest possible value of the induced 
number in this sector is — n/2. Any general state of the system is built from one of these 
states by additional integer number of fermions. All the states in the system therefore 
possess half-integral values for the fermion number if n is odd. 

One puzzle immediately arises, what is the meaning of negative values for the fermion 
number operator for Majorana fermions? In the trivial vacuum, we can identify the 
Majorana basis as 

V = + (6.13) 

This leads to the Majorana condition which results in identification of particles with 
anti-particles according to 

C^C f = ip (6.14) 



72 



making negative values for the number meaningless. Here C is the charge conjugation 
operator. We shall first verify that in the zero-mode sector we must indeed assign negative 
values to the number operator. It is sufficient to treat the case of a single zero-mode, which 
generalizes easily to any larger number of zero-modes. The number operator possesses 
the properties 

[N, ip] = -ip and [N, ift] = ^ (6.15) 



CNC ] = N (6.16) 

Had it been the Dirac case, there should be a minus sign on the right hand side of eq. 
(16.16)1 . This is absent due to the Majorana condition. The fermion field operator for the 
lowest winding sector is now expanded as 

)+ ^2b Ks u Ks (x) + h.c.\ (6.17) 

L k,s k,s ) 

where the first summation is over all the possible bound states of non-zero frequency with 
real space-dependence of the form ~ e - K X J- j n the transverse space directions xj_, and 
the second summation is over all unbound states, which are asymptotically plane waves. 
These summations are suggestive and their exact connection to the Weyl basis mode 
functions |102j are not essential for the present purpose. Note however that no "/i.e." is 
needed for the zero energy mode which is self-conjugate. Then the Majorana condition 
(16.14)1 requires that we demand 

C cC ] = c and C c f C ] = c f (6.18) 

Unlike the Dirac case, the c and are not exchanged under charge conjugation. The only 
non-trivial irreducible realization of this algebra is to require the existence of a doubly 
degenerate ground state with states |— ) and |+) satisfying 

c\-) = |+) and c f |+) = |-) (6.19) 

with the simplest choice of phases. Now we find 

C cC ] C\-) = C\+) (6.20) 
=► c(C|-» = (C|+» (6.21) 
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This relation has the simplest non-trivial solution 

C\-) = tjc\-) and C\+) = r]c\+) ( 6 - 22 ) 
where, for the consistency of (jSHSj) and ()6.21|) 77^ and 77J must satisfy 

(ifo)" 1 ^ = 1 (6.23) 

Finally we verify that we indeed get values ±1/2 for N. The standard fermion number 
operator which in the Weyl basis is 

N F = ^[*{* L -V L tf L ] (6.24) 

acting on these two states gives, 

~(cct - c t c ) |±) = ±I|±) (6.25) 

The number operator indeed lifts the degeneracy of the two states. For s number of 
zero modes, the ground state becomes 2 s -fold degenerate, and the fermion number takes 
values in integer steps ranging from — s/2 to +s/2. For s odd the values are therefore 
half-integral. Although uncanny, these conclusions accord with some known facts. They 
can be understood as spontaneous symmetry breaking for fermions |l(J3j . The negative 
values of the number thus implied occur only in the zero-energy sector and do not continue 
indefinitely to — oo. Instead of an unfathomable Dime sea we have a small Majorana pond 
at the threshold. 



6.4 Quantum mechanical stability 

The theory of eq. (J6.1j) possesses a gauge symmetry which is reflected in the effective the- 
ory (j6.4j) as Nl — > e ta NL, Nr — >e ia Npt giving rise to the usual conserved number for Dirac 
fermions. The lowest winding vortex sector results in half-integer values for this number. 
Quantum Mechanical stability of this sector follows from well known arguments HU4j 
which can now be understood as either following from distinctness of sectors of different 
values of (—1)^, or as a consequence of a residual subgroup of the gauge symmetry. For 
the Majorana case we shall now carry out this kind of argument explicitly. 
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It is known that Majorana fermions can be assigned a unique parity |l(J4j . either of 
the values ±i. Accordingly let us choose i to be the parity of the free single fermion states 
in the trivial vacuum. 

As a step towards deriving our superselection rule, we determine the parities of the 
zero-energy states. The fermion spectrum should look the same as trivial vacuum far 
away from the vortex. In turn the parities of the latter states should be taken to be the 
same as those of the trivial ground state. Next, any of these asymptotic free fermions 
is capable of being absorbed by the vortex (see for instance |lU5j ). In the zero energy 
sector this absorption would cause a transition from | — 1/2) to 1 1/2) and cause a change 
in parity by i. Thus the level carrying fermion number +1/2 should be assigned a parity 
e w/2 re ] a £i ve t the —1/2 state. Symmetry between the two states suggests that we assign 
parity e i7r/4 to the Np = 1/2 and e" i7r/4 to the N F = -1/2 states. 

Similar reasoning applies to a residual discrete symmetry belonging to the original U(l) 
gauge group of Lagrangian 1)6.11)1 . According to eq. ()6.13jl . under gauge transformation, 

^^a]^(e'^L + e-^) (6.26) 

Thus a = it preserves the choice of the Majorana basis upto a sign. After symmetry 
breakdown and Higgs mechanism, the Yukawa coupling takes the form ~ (m + 4>)ipip, 
which is invariant under the residual Z2 symmetry ip — > —ip. We can use this as a discrete 
symmetry distinguishing states of even and odd Majorana fermions. Since single Majorana 
fermions can be absorbed by the vortex )lU5j . the ground states |±) are distinguished from 
each other by a relative negative sign. To be symmetric we can assign the value ±2 to 
these states under this discrete symmetry with sign same as in the value of the number 
operator. It is possible to prove the superselection rule using this conserved quantity. 
However we also see that this discrete symmetry can be used to change our convention of 
the parity for free Majorana particles from +i to —i. Thus the two are intimately related 
and in what follows we shall use the parity with convention as in the preceding paragraph. 

We now show the inappropriateness of superposing states of half-integer valued fermion 
number and integer valued fermion number jH2]- The operation V , parity transformation 
performed four times must return the system to the original state, upto a phase. Consider 
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forming the state ^5 =^=(|l/2) + |1)) from states of half-integer and integer value for the 
fermion number. But 

VUs = -^(-|l/2) + |l» (6.27) 

Thus this operation identifies a state with another orthogonal to it. Similarly, application 
of V which should also leave the physical content of a state unchanged results in yet 
another linearly independent state, ^=(i|l/2) — |1)). Thus the space of superposed states 
collapses to a trivial vector space. The conclusion therefore is that it is not possible to 
superpose such sectors. In turn there can be no meaningful operator possessing non-trivial 
matrix elements between the two spaces. This completes our proof of the theorem. 
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Chapter 7 

B-L cosmic strings and baryogenesis 
via leptogenesis 

7.1 Introduction 

Cosmic strings, which are one dimensional topological defects, are expected to be formed 
during the early Universe phase transitions. When a larger symmetry group G breaks 
to H and if the corresponding vacuum manifold Ai = G/H is not simply connected (for 
definition see Appendix C) then the formation of cosmic strings is assured. 

The presence of cosmic string in the real physical space can be ensured by encircling it 
with a closed path (for definition see Appendix C) or loop f(t) as shown in the figure 17.11 
The field values 4>(x, t) at points along the loop, /(£), take values in Ai, and so they map 
points from a physical space R 3 into the vacuum manifold Ai. Hence <f>(x,t) completes a 
mapping from S 1 onto a path g(t) in Ai, g(t) = <p (f(t)). If the path g(t) possesses a non- 
trivial winding n in Ai, then a string must be present in the physical space. But this is 
exactly the criterion for g(t) to belong to the homotopy class n. This correspondence is one 
to one since each type of string can be identified with a unique element of vr 1 (S' 1 ) = Z. One 
might suppose, therefore, that string solutions in general can be classified by the elements 
of it\(Ai). Thus in general if TTi(Ai) 7^ / then the Higgs field <fi{x) will necessarily wind 
around Ai in a nontrivial way and a network of cosmic strings is formed. 

If G is a local gauge group then the strings formed after the breaking of G harbor 
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Figure 7.1: A non-zero winding in the phase ensures a cosmic string within. 



quanta of massive gauge bosons as well as zero-modes of right handed Majorana fermions 
[HTij . Usually the same Higgs which forms the string also gives mass to the Majorana 
fermions. Presence of right handed Majorana zero modes on the string leads to a curious 
phenomenon of fermion number fractionalization. However, the existence of the number 
of zero modes on a string depends on its winding number [ST]. If the winding number of 
the string is odd then the induced fermion number on the string is half integral. Hence the 
decay of a cosmic string to fundamental particles is forbidden because of the unavailability 
of quanta of particles carrying fractional fermion number in the translational invariant 
vacuum. In this chapter we shall study the decay of strings by forming loops since the 
string loops may not carry any fractional fermion number. 

In the context of type-I models, we study the formation, evolution and the decay of B — 
L cosmic strings. When the string loop decays it emits those harbored particles, which are 
collectively called X-particles. The delayed decay of these massive X-particles provides a 
link between the early Universe and recent cosmology. In particular, we consider the case 
of baryogenesis via the route of leptogenesis. Assuming the hierarchical mass spectrum 
in the right handed neutrino sector we discuss the leptogenesis constraint on the mass 
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scale of lightest right handed Majorana neutrino as well as the scale of B — L symmetry 
breaking phase transition. 

7.2 C/(1)b-l cosmic strings and neutrino zero modes 

An inevitable feature of U(l) gauge symmetry breaking phase transition is the formation of 
cosmic strings. In the vacuum manifold, M. = (77(1)//), since it\(M.) ^ I, a loop can not 
be trivially shrunk to a point and thus the existence of strings is assured. There are several 
realistic particle physics models where a gauged [/(1)b-l symmetry exists and breaks 
at a certain scale. Since SO(10) minimally incorporates [/(1)b-l gauge symmetry we 
consider the models embedded in 50(10). The following breaking schemes can potentially 
accommodate cosmic strings. One of the breaking schemes, motivated by supersymmetric 
50(10) [10(H 1107] . involves the intermediate left-right symmetric model: 

50(10) 54 + 45 SU{3) C ® SU(2) L ® SU(2) R ® C/(1) B -l 
126 + 126 SU{3) c ®SU{2) L ®U{l) Y ®Z 2 
10 + 10' SU(3) C ®U(1) Q ®Z 2 . (7.1) 

During the first phase of symmetry breaking, presumably at a GUT scale of ~ 10 16 GeV, 
monopoles are formed. However, during the second and third phases of symmetry breaking 
cosmic strings are formed since tti(^^) = 7r i(^f) = Z 2 , where the numbers inside the 
parentheses symbolize the group structures. The monopole problem in this model can 
be solved by using a hybrid inflation ending at the left-right symmetric phase of the 
Universe |107j thus inflating away the monopoles. The formation of cosmic strings in the 
later phases is of great interest since these "light" (i.e., lighter than GUT scale) cosmic 
strings do not conflict with any cosmological observations. The Z 2 strings of the low 
energy theory was investigated in an earlier work [80J. 

Another scheme is to break supersymmetric 50(10) directly to SU(3)c <S> SU{2)l ® 
^(1).r ® U(1)b-l with the inclusion of extra 54 + 45. The first 45 acquires a vacuum 
expectation value along the direction of B — L. However, the latter 45 acquires a vacuum 
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expectation value along the direction of T 3R : 



SO(10) 54 + 45 + 54' + 45' SU(3) C ® SU{2) L ® U{l) R ® £7(1) 



126 + 126 SU(3) C ® SU(2) L ® U(1) y ® Z 2 
10 + 10,' SU(3) C ® £7(1) Q ® Z 2 



(7.2) 



For the demonstration purpose it is sufficient to consider a model based on the gauge 
group SM eg) £/(1)b-l which is spontaneously broken to S'M. Existence of cosmic strings 
and the related neutrino zero-modes in this model can be established as follows. Let the 
gauge field corresponding to the £7(1) B _ L symmetry be denoted by C^, and the symmetry 
be broken by a SM singlet x- Let (x) be t^b-l below the critical temperature Tb-l- In 
a suitable gauge a long cosmic string oriented along the z-axis can be represented (in 
cylindrical polar coordinates) by the ansatz jl^l 



where n is an integer giving the winding number of the phase of the complex Higgs field 
X, and a is the gauge coupling constant for the group £7(1) B _ L . In order for the solution 
to be regular at the origin we set /(0) = g(0) =0. Also requiring the finiteness of energy 
of the solution, we set /(r) = g(r) = 1 as r — > oo. It turns out that both f(r) and 
g(r) take their asymptotic values everywhere outside a small region of the order of Tj^ L 
around the string. Thus away from the string (x) = ^b-l up to a phase 9, and C M is 
a pure gauge. The arbitrary phase 9 = 9(x) can vary in different regions of space. For 
X to be single valued 9 must change by an integer multiple of 27r around a closed loop. 
When the loop is shrunk to a point, 9 becomes undefined, so that there exists point 
where (x) = and the symmetry is unbroken i.e. thin tubes of false vacuum get trapped 
some where inside the loop. It is in this sense that the object is referred to as a 'defect', 
a region of unbroken symmetry (false vacuum) surrounded by broken symmetry region 
(true vacuum). Such strings are either infinitely long or closed loops. While an infinite 
straight string is topologically stable due to the non trivial winding of the Higgs, a closed 
string loop can decay by emitting gravitational radiation. The strings trap excess energy 



C, 




(7.3) 



(7.4) 
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density associated with gauge and Higgs fields, and, in addition, fermion zero modes which 
make the strings massive. The mass scale of the string is fixed by the energy scale of the 
symmetry breaking phase transition t/b-l at which the strings are formed. Then the mass 
per unit length of a cosmic string, fj,, is of order ?7 B _ L ~ Tg_ L . 
The Lagrangian for the right-handed neutrino is 

Cu* = m^D^VR - ~ [ifVUxvR + H.C\ , (7.5) 

where h is the Yukawa coupling constant, cr M = (— J, cr l ), and v c R = io 2 v R defines the Dirac 
charge conjugation operation. The resulting equations of motion have been shown (50] to 
possess |n| normalisable zero- modes in winding number sector n. The field equations in 
the U{1) example are 



-e»[d r + & + d z + d, 

d z -d t e- ld [d r -id ti 2i 



id r a i a ng(r) 



v R -M R e ine iy* R = 0, (7.6) 



where the expressions (17. 3 J) and ()7.4)1 have been substituted for x an d C^, and Mr = 
/^b-l- In the winding number sector n the normalizable zero-modes obey <j z vr = vr and 
are of the form 

Mr, 0) = ^ j (U(r)e iW + V*^^- 1 ^ 9 ) , (7.7) 

where U(r) and V(r) are well behaved functions at the origin and have the asymptotic 
behavior ~ exp(— Mrt)/^. When nontrivial z and t dependences are included, these 
modes have solutions that depend on z + t and are Right movers. For n < 0, normalizable 
solutions obey <j 3 i>r = —ur, and form the zero-energy set of a Left moving spectrum. On 
a straight string these modes are massless. However on wiggly strings they are expected 
to acquire effective masses proportional to the inverse radius of the string curvature. 
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7.3 Evolution of cosmic strings: Formation and evo- 
lution of closed loops and production of massive 
particles 

7.3.1 Scaling solution and closed loop formation 

The evolution of cosmic strings in the expanding Universe has been studied extensively, 
both analytically as well as numerically; for a text-book review, see the monograph [46 . 
Here we briefly summarize only those aspects of cosmic string evolution that are relevant 
for the present purpose, namely the formation and subsequent evolution of closed loops of 
strings and production of massive particles from them. This closely follows the discussion 
in section 6.4 of Ref. [TU3j . 

Immediately after their formation at a phase transition, the strings would in general be 
in a random tangled configuration. One can characterize the string configuration in terms 
of a coarse-grained length scale £ s such that the overall string energy density p s is given 
by p s = Initially, the strings move under a strong damping force due to friction with 

the ambient thermal plasma. In the friction dominated epoch a curved string segment of 
radius of curvature r acquires a terminal velocity oc 1/r. As a result the strings tend to 
straighten out so that the total length of the strings decreases. Thus the overall energy 
density in the form of strings decreases as the Universe expands. This in turn means that 
the length scale £ s increases. Eventually, l; s reaches the causal horizon scale ~ t. After 
the damping regime ends (when the background plasma density falls to a sufficiently 
low level as the Universe expands), the strings start to move relativistically. However, 
causality prevents the length scale £ s from exceeding the horizon size ~ t. Analytical 
studies supported by extensive numerical simulations show that the subsequent evolution 
of the system is such that the string configuration reaches a "scaling regime" in which the 
ratio y = x remains a constant. Numerical simulations generally find the number x to 
lie approximately in the range ~ 0.4-0.7. This is called the scaling regime because then 
the energy density in the form of strings scales as, and remains a constant fraction of, 
the energy density of radiation in the radiation dominated epoch or the energy density of 
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matter in the matter dominated epoch both of which scale as t~ 2 . 

The fundamental physical process that maintains the string network in the scaling 
configuration is the formation of closed loops which are pinched off from the network 
whenever a string segment curves over into a loop, intersecting itself. It is difficult to 
calculate ab initio the length distribution of the closed loops so formed, but numerical 
simulations find that loops are formed typically on the scale of the smallest scale structure 
allowed by the spatial resolution of the simulation. It is expected ^H] that the smallest 
scale structure on the string at any time time t would be ~ TGfj,t, which is determined by 
gravitational radiation from the small scale structure and its back-reaction on the string. 
Here T ~ 100 is a geometric factor [IB]- Thus, the closed loops at birth can be assumed 
to have a typical length j3H! 



and they are formed at a rate (per unit volume per unit time) which, in the radiation 
dominated epoch, is given by 



where T ~ 100 is a geometrical factor that determines the average loop length, and K is 
a numerical factor of order unity. 

The whole string network consisting of closed loops as well as long strands of strings 
stretched across the horizon gives rise to density fluctuations in the early Universe which 
could potentially contribute to the process of formation of structures in the Universe. More 
importantly, they would produce specific anisotropy signatures in the cosmic microwave 
background (CMB). Using a large-scale cosmic string network simulation and comparing 
the resulting prediction of CMB anisotropics with observations, a recent analysis P-10J puts 
an upper limit on the fundamental cosmic string parameter /i, giving Gfi <, 0.7 x 10 -6 . 
This translates to an upper limit, ffe-L ~ 1-0 x 10 16 GeV, on the symmetry-breaking 
energy scale of the cosmic string-forming phase transition. This probably rules out cosmic 
string formation at a typical GUT scale ~ 10 16 GeV. However, lighter cosmic strings 
arising from symmetry breaking at lower scales, such as the B — L cosmic strings in the 
case of the .SO (10) model discussed in the previous section, are not ruled out by CMB 



L b = KTGfit , 




drib 
~dJ 



(7.9) 
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anisotropy constraints. For detailed discussions of constraints on topological defects in 
general and cosmic strings in particular from CMB anisotropies see, for example [111J. 

It should be noted here that, in the standard scenario of cosmic string evolution 
described above, the loops are formed on a length scale that is a constant fraction of the 
horizon length, as given by equation (|7.8|) . Thus, the average size of the newly formed 
loops increases with time. At the relevant times of interest, these loops, although small 
in comparison to the horizon scale, would still be of macroscopic size in the sense that 
they are much larger than the microscopic string width scale w ~ Vb-l ~ ^ l ^ 2 - 

In contrast, results of certain Abelian Higgs (AH) model simulations of cosmic string 
evolution |112j seem to indicate that scaling configuration of the string network is main- 
tained primarily by loops formed at the smallest fixed length scale in the problem, namely, 
on the scale of the width w ~ Vb-l ~ ^ l ^ 2 of the string. These microscopic "loops" 
quickly decay into massive particles (quanta of gauge bosons, Higgs bosons, heavy fermions 
etc.) that "constitute" the string. In other words, in this scenario, there is essentially no 
macroscopic loop formation at all; instead, the scaling of the string network is maintained 
essentially by massive particle radiation. In order for the scaling configuration of the 
string network to be maintained by this process, the microscopic loops must be formed 
at a rate 



The above scenario of cosmic string evolution in which massive particle radiation rather 
than gravitational radiation plays the dominant role is, however, currently a subject of 
debate [113J. One of the major problems hindering a resolution of the issues involved is 
the insufficient dynamic range possible in the currently available AH model simulations 
and the consequent need for extrapolation of the simulation results to the relevant cos- 
mological scales, which is not straightforward. In the present case, we shall primarily 
restrict ourselves to consideration of the "standard" macroscopic loop formation scenario 
described by equations ()7.8|) and ()7.9|) above, although we shall have occasions to refer to 
the massive particle radiation scenario below (see, in particular, section 7.3. 2. B). 




(7.10) 
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7.3.2 Fate of the closed loops and massive particle production 

The behavior of the closed loops after their formation may be broadly categorized into 
following two classes: 

A. Slow death 

Any closed loop of length L in its center of momentum frame has an oscillation period 
L/2 j!14j . However, a loop may be either in a self- intersecting or non-selfintersecting con- 
figuration. In general, a closed loop configuration can be represented as a superposition 
of waves consisting of various harmonics of sin's and cos's. Some explicit low harmonic 
number analytical solutions of the equations of motion of closed loops representing non- 
selfintersecting loops are known in literature |114t 11151 11161 IH7| , and it is possible that 
there exists a large class of such non-selfintersecting solutions. Indeed, numerical simu- 
lations, while limited by spatial resolution, do seem to indicate that a large fraction of 
closed loops are born in non-selfintersecting configurations. 

A non-selfintersecting loop oscillates freely. As it oscillates, it loses energy by emitting 
gravitational radiation, and thereby shrinks. When the radius of the loop becomes of the 
order of its width w ~ Vb-l ~ A* 1 ^ 2 ) the loop decays into massive particles. Among these 
particles will be the massive gauge bosons, Higgs bosons, and in the case of the B — L 
strings, massive right-handed neutrinos (ur) which were trapped in the string as fermion 
zero modes. We shall hereafter collectively refer to all these particles as X particles. We 
are, of course, interested here only in the z/r's. In addition to those directly released from 
the loop's final decay, there will also be some z/r's coming from the decays of the gauge 
and Higgs bosons released in the final loop decay. It is difficult to calculate exactly the 
total number of ur's so obtained from each loop, but we may expect that it would be a 
number of order unity. For the present purpose we shall assume that each final demise of 
a loop yields a number Nn ~ 0(1) of heavy right handed Majorana neutrinos; we shall 
keep this number Nn free parameter. 

The rate of release of vr's at any time t by the above process can be calculated as 
follows. The lifetime of a loop of length L due to energy loss through gravitational wave 
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radiation is 

TGw-iTGnr'L. (7.11) 

Equations ()7.8j) and (|7.11J) thus show that loops born at time t have a lifetime ~ Kt >, iJ _1 (t), 
where i? _1 (t) ~ t is the Hubble expansion time scale. It is thus a slow process. From the 
above, we see that the loops that are disappearing at any time t are the ones that were 
formed at the time (K + Taking into account the dilution of the number density of 

loops due to expansion of the Universe between the times of their birth and final demise, 
equation ()7.9j) gives the number of loops disappearing due to this "slow death" (SD) pro- 
cess per unit time per unit volume at any time t (in the radiation dominated epoch) as 

= •/ SD ^2 UW^ ^ * =/sd(A+1) / — , (7.12) 

where /sd is the fraction of newly born loops which die through the SD process. 

The rate of release of the heavy right-handed neutrinos (we shall hereafter denote it 
by iV; see section 7.4 below) due to SD process can then be written as 

In a comoving volume the above injection rate (|7.1Hjl of massive Majorana neutrinos is 
given by 

fdYff\ 1.57 xlfT 1 ' / M N \ 2 ( M N 

\-w) SD z; — NNfsD {^rJ [Cev, 

where Y$ = n^/s is the density of massive Majorana neutrino in a comoving volume with 
s = 43.86g*T 3 being the entropy density and Z = M^/T is the dimensionless variable 
with respect to which the evolution of the various quantities is studied. Here we have 
used the numerical values for the constants V = 100, x = 0.5, g* = 100 and K — 1. We 
shall use this equation while solving the Boltzmann equations numerically in section 7.5. 
For our numerical purpose we shall be interested in the lightest right handed neutrino Ni 
whose mass is given by M x . Since Mi relates to the scale of symmetry breaking phase 
transition as Mi = /i^b-l hence for all values of fi < 1, M 1 in general satisfies the 
constraint M 1 < r\ B _ L . 
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B. Quick death 

Some fraction of the loops may be born in configurations with waves of high harmonic 
number. Such string loops have been shown |118j to have a high probability of self- 
intersecting. Ref. |118j gives the self-intersecting probability of a loop as 



si 



where a = 0.4, (3 = 0.2, and n is the harmonics number. 

A self-intersecting loop would break up into two or more smaller loops. The process 
of self-intersection leaves behind "kinks" on the loops, which themselves represent high 
harmonic configurations. So, the daughter loops would also further split into smaller 
loops. If a loop does self- intersect, it must do so within its one oscillation period, since 
the motion of a loop is periodic. Under this circumstance, since smaller loops have smaller 
oscillation periods, it can be seen that a single initially large loop of length L can break 
up into a debris of tiny loops of size %_ L (at which point they turn into the constituent 
massive particles) on a time-scale ~ L. Equation ()7.8|) then implies that a loop born at 
the time t in a high harmonic configuration decays, due to repeated self-intersection, into 
massive particles on a time scale tqd ~ KTGfit <C It is thus a "quick death" 

(QD) process — the loops die essentially instantaneously (compared to cosmological time 
scale) as soon as they are formed. Equation ()7.9|h therefore, directly gives the rate at 
which loops die through this quick death process: 

dn QD dn b 
dt = ' (7 - 16) 

where /qd is the fraction of newly born loops that undergo QD. 

Note that, since these loops at each stage self-intersect and break up into smaller loops 
before completing one oscillation, they would lose only a negligible amount of energy in 
gravitational radiation. Thus, almost the entire original energy of these loops would 
eventually come out in the form of massive particles. 

Assuming again, as we did in the SD case, that each segment of length ~ w ~ fi^ 1 ^ 2 
of the loop yields a number ~ 0(1) of heavy right-handed Majorana neutrinos, we 
can write, using equations (|7.16j) . (|7.9|) and (|7.8|) . the rate of release of the iV's due to 
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QD process as 

(If)^^ 172 '- 3 - (7 ' 17) 

In a comoving volume the above rate (J7.17|) can be rewritten as 

U^J QB " is ^MavKavj' (7 ' 18) 

We shall use this equation while solving Boltzmann equations numerically in section 7.5. 
However, for the leptogenesis purpose we shall be interested in the lightest right handed 
neutrino Ni and hence we replace by N\ and correspondingly Z = Mi/T. 

It is interesting to note here that if all loops were to die through this QD process, 
i.e., if we take /qd = 1 in equations (j7.16|) and (|7.17|) . then the situation is in effect 
exactly equivalent to the microscopic loop formation scenario described by equation (|7.10jl . 
although the primary loops themselves are formed with macroscopic size given by equation 

While the important issue of whether or not massive particle radiation plays a dom- 
inant role in cosmic string evolution remains to be settled, the standard model may, of 
course, still allow a small but finite fraction, /qd -C 1, of quickly dying loops. There 
already exist, however, rather stringent astrophysical constraints |119[ HU9j on /qd from 
the observed flux of ultrahigh energy cosmic rays (UHECR) above 10 11 GeV [120J and the 
cosmic diffuse gamma ray background in the energy region 10 MeV - 100 GeV measured 
by the EGRET experiment |121j . This comes about in the following way: 

The massive X particles released from the string loops would decay to SM quarks 
and leptons. The hadronization of the quarks gives rise to nucleons and pions with 
energy up to ~ Mx, the mass of the relevant X particle. The neutral pions decay to 
photons. These extremely energetic nucleons and photons, after propagating through the 
cosmic radiation background, can survive as ultrahigh energy particles. The observed flux 
of UHECR, therefore, puts constraints on the rate of release of the massive X particles, 
thereby constraining /qd- The most stringent constraint on /qd, however, comes from the 
fact that the electromagnetic component (consisting of photons and electrons/positrons) 
of the total energy injected in the Universe from the decay of the X particles initiates an 
electromagnetic cascade process due to interaction of the high energy electrons/positrons 
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and photons with the photons of the various cosmic background radiation fields (such as 
the radio, the microwave and the infrared/optical backgrounds); see, e.g., Ref. [Tp9j for a 
review. As a result, a significant part of the total injected energy cascades down to lower 
energies. The measured flux of the cosmic gamma ray background in the 10 MeV - 100 
GeV energy region |121j then puts the constraint [109J 

/QD^e^O.exlO- 6 , (7.19) 

where 7716 = (^b-l/10 16 GeV). For GUT scale cosmic strings with tjiq = 1, for exam- 
ple, the above constraint implies that /qd < 10~ 5 , so that most loops should be in 
non-selfintersecting configurations, consistent with the standard scenario of cosmic string 
evolution. Note, however, that /qd is not constrained by the above considerations for 
cosmic strings formed at a scale 7/b-l 3.1 x 10 13 GeV. 

In this context, it is interesting to note that there is no equivalent constraint on 
the corresponding parameter /sd for the slow death case from gamma ray background 
consideration. The reason is that, unlike in the QD case where the entire initial energy 
of a large loop goes into X particles, only ~ one X particle is released from a initially 
large loop in the SD case. This in turn makes the time dependence of the rate of release 
of massive particles oc t~ 4 in the SD case (see equation (J7.13)) ). while it is oc t -3 in the 
QD case (see equation (j7.17|0 . Thus, while the SD process dominates at sufficiently early 
times, the QD process can dominate at relatively late times and can potentially contribute 
to the non-thermal gamma ray background. 

7.4 Analytical estimation of baryon asymmetry 

7.4.1 Decay of heavy right-handed Majorana neutrinos and L- 
asymmetry 

We follow the same convention as given in section 2.3.1. In this scenario the lepton asym- 
metry is produced by the decay of right-handed heavy Majorana neutrino to SM lepton 
(£) and Higgs (0) through the Yukawa coupling. We assume a normal mass hierarchy in 
the right handed heavy Majorana neutrino sector, Mi < M 2 < M 3 . In this scenario it is 
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reasonable to expect that the final lepton asymmetry is produced mainly by the decay of 
the lightest right handed neutrino AV As the Universe expands, the temperature of the 
thermal plasma falls. Below a temperature Tp ~ M\, all L- violating scatterings mediated 
by Ni freeze out, thus providing the out-of-equilibrium situation necessary for the sur- 
vival of any net L-asymmetry generated by the decay of the AVs. The final L-asymmetry 
is, therefore, given essentially by the CP asymmetry parameter 

3 Mim 3 

ei = —° > ( 7 - 2 °) 

l07T IT 

where v ~ 174 GeV is the electroweak symmetry breaking scale and 8 < 1. 

An accurate calculation of the net L-asymmetry can only be done by numerically 
solving the full Boltzmann equations that include all lepton number violating interactions 
involving all the AVs present at any time, including the AVs of non-thermal origin such 
as the ones produced from the decaying cosmic string loops, as well as those of thermal 
origin. This will follow up our analytical estimations. For analytical calculation, first 
we shall simply assume that below the temperature Tp = Mi, all interactions except the 
decay of the Ni are unimportant, so that each Ni released from cosmic strings additively 
produces a net L-asymmetry when it decays. 

For the present purpose, we note that there is an upper bound [2Z1I2E] on e i 5 which is 
related to the properties of the light neutrino masses. In a standard hierarchical neutrino 
mass scenario with m 3 ^> m-i > mi, this upper limit is given by fZT\ l2*%] 



N = 1g _ o • (7-21) 



3 Mim 3 

167T V 2 

The above upper limit is in fact saturated j2H| in most of the reasonable neutrino mass 
models, which we shall assume to be the case. 

Assuming the standard light neutrino mass hierarchy, the heaviest light neutrino is 

1 /2 

given by ~ (Am 2 tm ) ~ 0.05 eV. In our calculations below, we shall use 

ei * 9.86 x 10- ( ^ - ) ( (Am H! /2 ' ] • (7.22) 



^10 13 GeVy I 0.05 eV 

The L-asymmetry is partially converted to a 5-asymmetry by the rapid nonperturba- 
tive sphaleron transitions which violate B+L but preserve B—L. Assuming that sphaleron 
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transitions are ineffective at temperatures below the electroweak transition temperature 
(Tew)j the 5-asymmetry is related to L-asymmetry by the relation [4*lj 

B = p (B — L) = — — L ~ — 0.55L , (7.23) 
p — 1 

where we have taken p = 28/79 appropriate for the particle content in SM |44j . 
The net baryon asymmetry of the Universe is defined as 

Y B = ?SZ2*, (7.24) 
s 

where s = 43.86g*T 3 is the entropy density, with being the number of relativistic 
degrees of freedom contributing to the entropy at the temperature T. At temperatures in 
the early Universe relevant for the process of baryon asymmetry generation, ~ 100 in 
SM. 

Observationally, the BAU is often expressed in terms of the baryon-to-photon ratio 
rj = (rie — ng)/n 7 , whose present-day- value rjo is related to that of Yb through the relation 

7/o~7.0y fliO . (7.25) 

We now proceed to estimate the contribution to the BAU from the two cosmic string 
loop processes discussed in the previous section. 

7.4.2 Slow death case 

The contribution of the SD process to go can be written as 

r? SD ~ 7.0 x 0.55 e, jT ^ dt , (7.26) 

where tp is the cosmic time corresponding to the temperature Tp ~ Mi and to is the 
present age of the Universe. Using equations (|7.13j) and the standard time-temperature 
relation in the early Universe, 

where Mpi ~ 1.22 x 10 19 GeV is the Planck mass, we see that the dominant contribution 
to the integral in equation (|7.26|) comes from the time tp t , i.e., from the epoch of 
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temperature T F ~ Mi, giving 

"» SD = 2 ^ 10 - 7 *(lo4v) 4 (lo^v)^ 

= 2.0 x 10-AV/f (i^gv) 2 . (7.28) 

where we have defined the Yukawa coupling /i = Mi/?7b-l, used equation ()7.22|) for ei 
with (Am^ tm ) 1//2 = 0.05 eV, and also taken x = 0.5, T = 100, if = 1 and /sd = 1 in 
equation (j7.13|) . 

The Yukawa couplings are generally thought to be less than unity. With fx < 1, we 
see from (|7.28j) and (J2.20)) that the cosmic string loop slow death process can produce the 
observed BAU only for B — L phase transition scale 

/7b-l ~ 5.5 x 10 n iV~ 1/2 GeV . (7.29) 

Assuming Nn <, 10, say, we see that cosmic string loop SD process can contribute to 
BAU if r]B-L ~ 1-7 x 10 11 GeV; lower values of 7/b-l are relevant only if we allow f\ > 1. 

At the same time, for a given //b-l satisfying ()7.29j) . in order that the contribution 
()7.28|) not exceed the highest allowed observed value of 770 given by equation ()2.20|) . the 
Yukawa coupling fi must satisfy the constraint 

which, in terms of the lightest heavy right handed Majorana neutrino mass Mi, reads 

M*> < 2.4 x 10 12 iV- 1/4 (j^y) V2 GeV . (7.31) 

Note the rather weak dependence of the above constraints on Nn- Also, the 4th power 
dependence on Mi of equation ()7.28|) and the rather narrow range of the observed value 
of 770 given by equation ()2.20|) together imply that, in order for the SD process to explain 
the observed BAU, Mi (and equivalently fi) cannot be much smaller than their respective 
values saturating the above constraints. 

It is important to note that in deriving the above constraint on Mi we assume that 
there is a strong hierarchy in the low energy neutrino sector. However, this limit becomes 
worse in the case degenerate neutrinos |122j 
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7.4.3 Quick death case 

Replacing (^f-) SD in equation ()7.26|) by (^f-) QD given by equation ()7.17|) . and following 
the same steps as in the SD case above, we get the contribution of the QD process to 770 
as 

* 5,7 xl 0-^ /QD (_4_) 2 (-^_) 

= 5,7 x IO^AV/qd/? (j^y)' ■ (732) 

From (J7.32j) and (|2.2()j) we see that, considering the most optimistic situation with /qd = 
1, the QD process is relevant for BAU only for 

ryg° L > 1.1 x lO 14 ^ 173 GeV; (7.33) 

lower values of 7/b-l are relevant only if we allow fx > 1. On the other hand, the constraint 
()7.19|) allows /qd = 1 on/y if ^b-l < 3.1 x 10 13 GeV. This can be reconciled with the 
above constraint ()7.33|) only for Nn > 45 or so. Such a large value of Nn seems unlikely. 
In general, using the constraint (|7.19jl on /qd in ()7.32j) we get 

Comparing again with the observed value of 7/0, we now see that, for /1 < 1, the QD 
process can be relevant for BAU only for 

rfH L > 1.2 x lO 15 ^ 1 GeV . (7.35) 

For values of t]b-l satisfying the above constraint (|7.35|) . the QD process can produce 
the observed value of BAU for 

which in terms of Mi now reads 

1 /2 

M X QD < 3.6 x lO 15 ^ 2 ( 1Q i 6 B G L eV ) GeV . (7.37) 

From the above discussions we see that, as far as their contributions to the BAU 
is concerned, the QD process becomes important only at relatively higher values of the 
symmetry breaking scale 7?b-l compared to the SD process. 
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7.5 Numerical calculation of baryon asymmetry 



In section 7.4, neglecting the washout effects we saw that the final lepton asymmetry, 
produced by the string emitted right handed neutrinos, depends on the density of Ni's 
at the freeze out epoch, Tp ~ M 1; as well as the amount of CP- violation due to the 
decays of Ni. Note that, the final lepton asymmetry does not depend on the source of 
iVi's whether they are generated by any thermal or non-thermal processes. Therefore, we 
attempt to establish a complete analysis of baryon asymmetry by including the iVi's of 
thermal origin as well of cosmic string origin . We then discuss the constraints on the 
scale of B — L symmetry breaking phase transition from the observed baryon asymmetry. 

7.5.1 Boltzmann Equations 

Considering the lightest right handed heavy Majorana neutrino (Ni) of cosmic string 
origin as well of thermal origin at any epoch of temperature T (or equivalently Z) the 
total rate of change of the abundance of iVYs is given by 

dY Nl ( dY Nl \ / dY Nl \ ^ 3g ^ 



dZ \ dZ J nQ V dZ / . . 

x ' D,b \ ' injection 



The first term on the right hand side of equation (|7.38|) is given by the usual Boltzmann 
equation ITT] 

= -(D + S){Y N1 -Y e N \) , (7.39) 



iVl 

dZ J D,S 

where D and S constitute the decay and AL = 1 lepton number violating scatterings 
which dilute the number density of N\, and Y^\ is the abundance of N\ in the thermal 
equilibrium situation. 

The second term on the right hand side of equation (J7.38)) gives the rate of injection 
of iVi's, coming from the disappearance of string loops, into a comoving volume at a rate 



+ ( d J?L\ j (7 _ 40) 



dZ I injection \ dZ J SD V dZ / QD 



where the two terms on the right hand side are given by the equations (j7.14j) and (j7.18j) 
respectively. One can see from the equations (|7.14j) and (j7.18|) that while the first term, 
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i.e, the slow death of string loops dominates at early times, the second term, i.e., the 
quick death of cosmic string loops dominates at late times. 

The two terms in equation ()7.38|) compete with each other. While the first term dilutes 
the density of iVi's in a comoving volume, the second term try to enrich it due to the 
continuous injection of N^s from the shrinkage of cosmic string loops. The dilution of the 
density of N^s, due to the CP- violating decay, vis-a-vis produces a net B — L asymmetry 
dynamically. This can be predicted by solving the Boltzmann equation 

HY 

= _ eijD (y m _ y«) - WY B _ L . (7.41) 

The first term in equation (|7.41j) . involving the decay term D, produces an asymmetry 
while a part of it gets erased by the wash out terms involved in W. Note that W 
includes the processes of inverse decay and the AL = 1, AL = 2 lepton number violating 
scatterings. The different terms D, S and W used in equations ()7.39j) and ()7.41|) are 
explained in section 2.3.3 and a summary of all the lepton violating scattering densities 
are given in section 4.3. Essentially the scattering densities depend on the two parameters 
Mi and rh\, the effective neutrino mass parameter. 

7.5.2 Constraint on effective neutrino mass (mi) 

In order to explain the present baryon asymmetry of the Universe all the lepton number 
violating interactions mediated by N\ have to satisfy the Sakharov's criteria pQ. While 
the complex nature of Yukawa coupling ()2.6|) violates CP at the same time the decay 
rate of Ni has to be less than the Hubble expansion parameter H in order to satisfy the 
out of equilibrium condition. This imposes an immediate constraint ()4.6|) on the effective 
neutrino mass parameter rh\ < 2 x 10 -3 eV as explained in section 2.3.3. Therefore, for all 
fhi < lO" 3 ^ at an epoch T < Mi a net B — L asymmetry can be generated dynamically. 

The above constraint on fhi can be envisaged in a model as follows. We assume a 
charge- neutral lepton symmetry |123j . In this scenario we take the texture of the Dirac 
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mass of the neutrino to be 



( n 'm e m n > 



V' 



nn> - v '///, ///,, ///,, ^in, (7.42) 
y y/m e m T m T . J 

Using ()7.42|) and ()2.22j) in equation ()2.28|) we get the constraint 

mi < 0.525 x 10~ 5 eV, (7.43) 
which is in concordance with our assumption rh\ < 1CT 3 eV at an epoch T < Mi. 

7.5.3 Solutions of Boltzmann equations 

At an epoch T ^> Mi the lepton number violating processes are sufficiently fast as to set 
the B — L asymmetry to zero. As the temperature falls and becomes comparable with Mi 
a net B — L asymmetry is generated through the CP- violating decay of N±. The resulting 
asymmetry can be obtained by solving the Boltzmann equations. We solve the equations 
(J7.38)) and (J7.41)) numerically with the following initial conditions 

Yft = Y°l znd Y™ L = 0. (7.44) 

Using the first initial condition we solve equation (|7.38|) for Y^, and the corresponding 
B — L asymmetry Yb-l is obtained from equation ()7.41|) by using the second initial 
condition of equation ()7.44|) . 

In the usual thermal scenario the B — L asymmetry depends on not only fhi, but also 
Mi. In the present case the B — L asymmetry depends, additionally, on 7/b-l since the 
injection rate of AVs from the cosmic string loops strictly depends on it. Depending on 
the B — L symmetry breaking scale, the effects of cosmic strings on the baryon asymmetry 
are shown in the Figures 17.21 and 1731 In these Figures we have taken the saturated value 
of the CP-asymmetry parameter e\. 

In Figure l?~2T a) it can be seen that in a thermal bath the B — L asymmetry approaches 
the final value at around Mi = 10T when all the wash out processes fall out of equilibrium. 
In contrast to it, in the presence of cosmic strings the B — L asymmetry continues to build 
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Figure 7.2: The evolution of B-L asymmetry with maximal allowed CP violation param- 
eter ei for different values of Mi, with rh\ = 10 _4 eV and i]b-l = 10 13 GeV (a) in absence 
of cosmic strings, and (b) in presence of cosmic strings. 




1000 




1000 



Figure 7.3: The evolution of B-L asymmetry with maximal allowed CP violation param- 
eter €\ for different values Mi, with rhi = 10~ 5 eV and t]b~l = 10 13 GeV (a) in absence of 
cosmic strings, and (b) in presence of cosmic strings. 
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up until the injection rate of AVs in a co moving volume is insignificant. For fi = l this 
happens around Mi = 6 x 10 2 T, as shown in Figure [T^b), which is far larger than in the 
purely thermal case. As a result of this, in presence of cosmic strings, for a fixed value of 
the symmetry breaking scale t/b-l = 10 13 GeV, the final B — L asymmetry is enhanced 
by three orders of magnitude for the effective neutrino mass fhi = 10~ 4 eV (Fig. I7.2lf b)) 
and by two orders of magnitude for fhi = 10~ 5 eV (Fig. 17.3( b)). For f\ < 0.01 the effect 
of cosmic string essentially disappears. 

The above happenings, in particular the dependence on the value of can be un- 
derstood as follows: In the absence of injection term when the wash out processes fall 
out of equilibrium the asymmetry produced by the decay of AVs does not get wiped out, 
and the produced B — L asymmetry remains as the final asymmetry. Since the decay 
rate of AVs depends linearly on fhi as inferred from equation ([4.10)1 . a relatively larger 
value of fhi implies that the condition of decay in out-of-equilibrium situation is satisfied 
at a relatively later time when the abundance of thermal AVs is relatively smaller, thus 
yielding a smaller final value of Yb-l- And as expected, this effect is larger for larger 
values of Mi . 

Now let us see what happens when the injection term due to the string loops is added. 
For values of T]b-l for which the string contribution remains subdominant to the thermal 
contribution, the dependence on fhi is essentially same as in the absence of the strings as 
explained above. However, the situation is reversed for those values of tjb-l for which the 
string contributions are dominant: Since the string-produced Aq's dominate over those of 
thermal origin at late times, the B — L asymmetry produced by the decay of the string- 
produced AVs automatically satisfy the out-of-equilibrium condition. In such a situation, 
a relatively larger value of fhi (implying a relatively larger rate of decay of the AVs) 
simply leads to a quicker rate of development of the Yb-l & t a relatively earlier time 
when the injection rate of the AVs is relatively larger, thus yielding a relatively higher 
final value of Yb-l- Figures IT^T b) and !7.3f b) plotted for fhi = 10~ 4 eV and fhi = 10~ 5 eV 
simultaneously bear out these expectations. 

Allowing a maximal CP-asymmetry, we see that for a reasonable value of fhi in certain 
models the string contribution overproduces the baryon asymmetry in comparison to the 
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thermal value for t]b-l = 10 13 GeV and fi >, 0.1. This suggests a lower 7/b-l or 5 < 1 
in equation ()7.20|) may be in better accord with the observed baryon asymmetry in the 
presence of B — L cosmic strings. Therefore, we have repeated the above calculations for 
t^b-l = 10 12 GeV and 7/b-l = 10 11 GeV. The results are summarized in tables ()7.5.3j) and 
(I7.5.3|) . In both cases one can see that for Yukawa coupling f\ < 1 and fh\ = 10 _5 eV the 
effect of cosmic strings disappear, i.e. (AY B _ L /Y^ h _ L ) — > 0, where AV B -l = Y-^X 8 * — ^b-l- 
This is in concordance with our assumption that rh\ < O(10 -5 ) eV in the given model in 
section 7.5.2. 

Studying the needed bound on S we find that for t/b-l = 10 13 GeV and Mi = 10 13 GeV 
to produce the final B — L asymmetry, Y^j^ 1 = O(10~ 10 ), we need 5 = O(10 -7 ) and 
O(10~ 6 ) for fhi = 10~ 4 eV and rh\ = 10~ 5 eV respectively. These values of 5 diminish the 
purely thermal contributions, Yb-l = O(10 -13 ) and O(10 -12 ) respectively. Thus there 
exists ranges of parameter values where, while the thermally abundant heavy neutrinos 
are not sufficient to produce the required B — L asymmetry, the cosmic string contribution 
can produce the required asymmetry. 



Table 7.1: Effect of cosmic strings on final B — L asymmetry is shown for different values 
of the Yukawa coupling f\ at t^b-l = 10 12 GeV 



fl 


m x = 10~ A eV 


m 1 = 10 5 eV 


Vth 
1 B-L 


-yth+st 
1 B-L 


log(5 reg(i ) 


Vth 
1 B-L 


\/-th+st 
1 B-L 


log(5 re(?rf ) 


1 


2.721 x 10~ 7 


1.231 x KT 5 


-5 


3.903 x 10~ 7 


1.594 x 10~ 6 


-4 


0.1 


3.098 x 10~ 8 


3.218 x 10~ 8 


-2 


3.970 x 10~ 8 


3.982 x 10" 8 


-2 


0.01 


3.254 x 10" 9 


3.254 x 10~ 9 


-1 


4.038 x 10" 9 


4.038 x 10~ 9 


-1 


0.001 


3.393 x 10" 10 


3.393 x 10~ 10 





4.107 x 10~ 10 


4.107 x 10- 10 






While the mechanism investigated here can be easily generalized to a supersymmetric 
model, there are unlikely to be significant changes in the quantitative aspects. In partic- 
ular unless some salient physics significantly modifies the scale of efficacy, viz., 10 11 GeV 
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Table 7.2: Effect of cosmic strings on final B — L asymmetry is shown for different values 
of the Yukawa coupling fi at ??b-l = 10 11 GeV 



h 


mj = 10" 4 e\/ 


m 1 = 10 5 eV 


Yb-l 


Yth-\-st 
B — L 


\og(5 reqd ) 


Yb-l 


Yth+st 
B — L 


log(<5 re<?d ) 


1 


3.098 x 10- 8 


1.514 x 10~ 7 


-3 


3.970 x 10- 8 


5.174 x -8 


-2 


0.1 


3.254 x 10~ 9 


3.266 x 10~ 9 


-1 


4.038 x 10~ 8 


4.039 x 10~ 9 


-1 


0.01 


3.392 x 10~ 10 


3.392 x 10~ 10 





4.107 x 10~ 10 


4.107 x 10~ 10 






of this mechanism, the supersymmetric generalization is not helpful in resolving the issue 
raised by the gravitino bound and the need for a mechanism at scales < 10 9 GeV. 

A more fruitful approach towards success of this scenario at lower energy scales would 
be to also seek additional sources for CP violation for the cosmic string generated heavy 
neutrinos. This is possible since they get produced from the decay of a bosonic condensate 
consisting of the string loops. However the strength of the present investigation is the 
direct comparison between the two sources of heavy neutrinos. If the mechanism appli- 
cable to the non-thermal source is different, the model becomes less constrained and its 
verifiability is sacrificed to a certain extent. 

A useful conclusion from the rather constrained scenario considered here is that in 
the regime ?7b-l > 10 11 GeV and hi >, 0.01 for ei <C e™ ax (i.e. <5< 1) while the thermal 
abundance of right handed neutrinos is not sufficient to produce the required B — L 
asymmetry the cosmic strings can give rise to the observed level of the asymmetry. 
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Chapter 8 



Conclusions 



Baryogenesis via leptogenesis is an attractive scenario that links the physics of right 
handed neutrino sector with the low energy neutrino data. In the light of current neu- 
trino oscillation data we studied the bounds on the mass scale of lightest right handed 
neutrino as well as their mass hierarchy from the leptogenesis constraint which we discuss 
in two different parts of the thesis. In part-I, we study the baryogenesis via leptogenesis 
in a thermal scenario, while part-II of the thesis is devoted to a study of the same in a 
non-thermal scenario. Moreover, we divide the phenomenological models into two cate- 
gories, type-I and type-II, depending on the seesaw mechanism used to generate the light 
Majorana neutrino masses. 

In Part-I, we begin with a brief introduction to type-I models. As an example, we 
consider a minimal extension of the SM with the gauge group SU(2) L x U(1) Y x U(1) Y ', 
where Y 1 is a linear combination of Y and B — L. Since B — L is a gauge symmetry of the 
model any pre-existing asymmetry is washed out. A net asymmetry is generated when 
B — L symmetry breaks. In such models, the right handed heavy Majorana neutrinos, iVj 
with % — 1, 2, 3, are singlet under the SM gauge group SU(2) L x U(1) Y ■ The type-I seesaw 
mechanism then gives rise to the light Majorana neutrino mass matrix, m v = m^ ) M^ i 1 m D) 
with mo is the Dirac mass matrix of the neutrinos. In these models, we briefly discussed 
the baryogenesis via the route of leptogenesis in a thermal scenario. The leptogenesis 
occurs via the out of equilibrium decay of thermally generated heavy Majorana neutrinos. 
We assume that the mass spectrum of right handed Majorana neutrinos is in normal 
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hierarchy. In this scenario, the final lepton asymmetry is given by the CP-violating 
decays of the lightest right handed neutrino Ni to SM Higgs and lepton. A part of this 
asymmetry is then transformed to the P-asymmetry through the equilibrated sphaleron 
processes. 

The L-asymmetry predicted by a model is a measure of the magnitude of CP- violation 
in that model. The theoretical upper bound on the CP-asymmetry, produced by the decay 
of A 7 !, in type-I models is given by e x < (3M 1 /16^^v 2 )^/Am^^, with Am^ im is the mass 
squared difference of light neutrinos in atmospheric data. The analytical estimation of 
L-asymmetry in these models recasts the upper bound on CP- asymmetry in terms of a 
lower bound on the mass scale of Ni to be Mi > 10 9 GeV. In the light of current neutrino 
oscillation data we numerically check the compatibility of the analytical ound on Mi by 
solving the required Boltzmann equations. 

The B — L gauge symmetry in the given example of type-I model is quite ad hoc. We 
therefore consider models where B — L gauge symmetry emerges naturally. In particular, 
we consider Left-Right symmetric model, a specific example of type-II seesaw models. 
The scalar sector of the model is very rich and consists of two triplets, namely Al and 
A R , and a bidoublet Higgs <3>, which contains two copies of SM Higgs. The type-II seesaw 
mechanism in this model gives rise the Majorana mass matrix of the light neutrinos of 
the form m v = Ml — m^M^m 1 ^- In contrast to type-I models, in the present case the 
additional mass, Ml, is provided by the vacuum expectation value of the triplet Al. The 
two terms, M L = m}J and m^M^m^ = ml, contributing to m u are called type-II and 
type-I respectively. In these models, irrespective of the magnitudes of type-I and type-II 
terms, we show that the lower bound on Mi can be reduced by an order of magnitude in 
comparison to type-I case. 

Within the Left-Right symmetric model, we consider a special case in which the CP- 
violation arises through the spontaneous symmetry breaking. The Lagrangian of the 
model is CP invariant which demands that all the Yukawa couplings should be real. In 
this scenario, the vacuum expectation values (VEVs) of the neutral Higgses are complex 
and they lead to complex masses for fermions and hence CP-violation. In the Left-Right 
symmetric model, there are four complex neutral scalars which acquire VEVs. However, 
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the unbroken global U(l) symmetries associated with SU(2) L and SU{2) R gauge groups 
allow two of the phases to be set to zero. Using the remnant U(l) symmetry after the 
breaking of SU(2) R , one phase choice is made to make the VEV of A R , and hence 
the mass matrix of right handed neutrinos, real. The phase associated with the other 
U(l) symmetry can be chosen to achieve two different types of simplification of neutrino 
mass matrix. In the type-II choice, the ml is made real leaving the CP- violating phase 
purely with m^ . In this phase convention, we derive a lower bound on the mass scale 
of iVi from the leptogenesis constraint by assuming a normal mass hierarchy in the right 
handed neutrino sector. It is shown that the mass scale of Ni satisfy the constraint 
Mi > 10 8 GeV, which is in good agreement with the lower bound on Mi in generic type-II 
seesaw models. In the type-I phase choice, only the type-I term contains CP-violating 
phase leaving type-II term real. This allows us to derive an upper bound on the heavy 
neutrino mass hierarchy from the leptogenesis constraint. In order to achieve the observed 
baryon asymmetry of the present Universe, it is found that the mass hierarchy of right 
handed neutrinos satisfy the constraint M 2 /Ml < 17 and M3/M1 < 289 simultaneously. 
Numerically we verified that these bounds are compatible with the low energy neutrino 
oscillation data for all values of Mi > 10 8 GeU as implied by the lower bound on M 1 
in type-II phase convention. The bound on Mi, in production of L-asymmetry through 
the CP- violating decays of thermally generated Aq, is > 10 8 GeU which is far above the 
current accelerator energy range and beyond the reach of the next generation accelerators. 
However, these scenarios are well motivated by the current status of low energy neutrino 
oscillation data. 

As an alternative, we consider mechanisms which work at TeV scale and be consistent 
with the low energy neutrino data. We study a general scenario for TeV scale leptogenesis 
in a gauged B — L symmetric model. By solving the Boltzmann equations we explore the 
viable regions in the plane of the effective neutrino mass arameter rh\ and the mass of 
lightest right handed neutrino Mi. We assume that the required lepton asymmetry of the 
present Universe is produced during the B — L gauge symmetry breaking phase transition. 
The limited erasure by the processes mediated by Aq requires that rh\ < lCT 3 eU. Solution 
of the relevant Boltzmann equations shows that for rh\ = !CT 4 eU the mass of Aq has to 
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be smaller than 10 12 GeV and can be as low as lOTeV. In a more restrictive scenario 
where the neutrino Dirac mass matrix is identified with that of the charged leptons it is 
necessary that Mi > 10 8 GeU in order to satisfy rh\ < !0~ 3 eV. In this scenario all values 
Mi, 10 8 GeU < Mi < 10 12 GeV can successfully create the required asymmetry. If the 
Dirac mass scale of neutrinos is less restricted, much lower values of Mi are allowed. In 
particular, a right handed neutrino as low as lOTeV is admissible if the Dirac mass scale 
is about a factor of 100 smaller than the mass scale of charged leptons. We conjecture 
that the hypothesis of TeV scale right handed neutrinos can be verified in the near future, 
thus providing an indirect evidence of baryon asymmetry creation at TeV scale. 

In part-II of the thesis, we discuss the formation and evolution of topological defects, in 
particular cosmic strings, in the early Universe. Topological defects arise as the solitonic 
solutions in gauge theories. 'Solitons' or 'solitary waves' are the time independent solu- 
tions of non-linear wave equations in classical field theories. The prime among them is A0 4 
theory. In 1+1-dimensions the solitonic solutions in A0 4 theory are called 'kinks'. Each 
solitonic solution is designated by a number called the 'topological charge' or 'winding 
number'. 

An inevitable feature of the early Universe phase transitions is the formation of topo- 
logical defects. In particular, we deal with cosmic strings. These defects are extended 
objects and are not distributed thermally. Therefore, the decay of these objects can be 
a non-thermal source of massive particles that constitutes them. Moreover, the cosmic 
strings formed at a phase transition can also influence the nature of a subsequent phase 
transition that may have important implications for the generation of Baryon Asymmetry 
of the Universe (BAU). 

In Quantum Field Theory (QFT), solutions of Dirac equation in the presence of 
solitonic objects lead to a curious phenomenon of 'fractional fermion number'. This is 
because of the existence of degenerate zero energy modes of fermions while quantized in 
the background of a solitonic vacuum. In contrast to it, in the translational invariant 
vacuum there are no zero energy solutions of Dirac equation and therefore fermions are 
quantized by integral unit. The fractional solitonic states are therefore superselected from 
the normal vacuum and are not allowed to decay in isolation. 
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An important feature of cosmic strings is that during their formation they trap zero 
energy modes of fermions. These fermionic zero modes induce fractional fermion number 
(|n|/2) on a string of winding number n. If n is odd then the induced fermion number 
on the string is half-integral. Therefore, it is superselected from the translation invariant 
vacuum where the eigenvalues of the number operator carry integral fermion number. 
Thus, a string of half integer fermion number can not decay in isolation because of the 
unavailability of fractional fermion states in the translational invariant vacuum. We con- 
struct examples where metastable infinitely long cosmic strings enjoy stability due to this 
quantum mechanical phenomenon. 

There exist both analytical as well as numerical studies of the evolution of cosmic 
strings network in the early Universe. These suggest that the strings network quickly 
enters a scaling regime in which the energy density of the strings scales as a fixed fraction 
of the energy density of radiation in the radiation dominated epoch, or the energy density 
of matter in the matter dominated epoch. In both cases the energy density scales as 
t~ 2 . In this regime one of the fundamental physical process that maintains the strings 
network to be in that configuration is the formation of sub-horizon size closed loops which 
are pinched off from the network whenever a string segment curves over into a loop, 
intersecting itself. 

There have been many scenarios that consider decaying, collapsing, or repeatedly 
self-intersecting closed loops of such cosmic strings providing a non-thermal source of 
massive particles that "constitute" the string. The decay of these massive particles gives 
rise to the observed B- asymmetry or a significant fraction of it. We have estimated the 
contribution to BAU from cosmic string loops which disappear through the process of 
(a) slow shrinkage due to energy loss through gravitational radiation — which we call 
slow death (SD), and (b) repeated self- intersections — which we call quick death (QD). 
We find that for reasonable values of the relevant parameters, the SD process dominates 
over the QD process as far as their contribution to BAU is concerned. We assume 
that the final demise of each string loop in the SD process produces 0(1) right handed 
neutrino, N\. We demand that the 5-asymmetry, produced by the decay of non-thermally 
generated iVi's through the leptogenesis route, should not exceed the observed value 
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predicted by WMAP. This requires that the mass scale of Ni satisfy the constraint Mi < 
2.4 x 10 12 (i] B ^ L /10 n GeV) 1/2 GeV, where i] B _ L is the scale of E/(1) b _l gauge symmetry 
breaking phase transition. This in turn constrains the B — L symmetry breaking scale to 
be 7/b-l > 5.6 x 10 11 GeV r unless the Majorana Yukawa coupling of Aq is allowed to be 
greater than unity. In this analytical approximation we have assumed that above the mass 
scale of N 1 there is no lepton asymmetry. A net lepton asymmetry has been produced 
just below the mass scale of Aq by it's CP-violating decays to SM Higgs and lepton. 
We then checked the analytical results against the numerical smulations. This is done by 
solving the relevant Boltzmann equations and including the effects of both thermal and 
string generated right handed neutrinos. We explored the parameter region spanned by 
the relevant light- and heavy neutrino mass parameters rh\ and Mi, and constrained the 
scale of B — L symmetry breaking, i] B _ L , as well as the Majorana Yukawa coupling f\ of 
the lightest right handed neutrino. It is shown that for the values t\b-l > 10 n GeV, where 
they can be effective, cosmic strings make a more dominant contribution than thermal 
leptogenesis. This, in turn, provides an upper bound on the CP violating phase giving 
rise to the L-asymmetry. 
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Appendix A 



The Boltzmann equations 



A.l Particle distribution in thermal plasma 

Let fi(p,t) be the density distribution of the particle species % in phase space. Assuming 
that the Universe is isotropic fi can be chosen to be independent of the position vector 
r. In the high temperature thermal bath it is assumed that all species of particles in 
the Universe were initially in thermal equilibrium and spread homogeneously. The only 
deviation comes for the massive particles as the Universe expands. The expansion of 
the Universe causes the momenta of all particles to redshift, so that p ~ So long as 
the energy density of the Universe is dominated by the ultra relativistic particle species 
the temperature of the Universe will like wise redshift as T ~ j^. Thus the distribution 
function, fi, is simply given by 

fi(p) = e-* (A.l) 

where % stands for the particle species. However, for massive particles the distribution 
function must be 

m=e- y -^r • (A.2) 

In the configuration space the number density of the of the particle species i is 

m = 9i J ( A - 3 ) 
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where g>j is the number of accessible spin states for the particle species i. Thus for a 
massless thermal photon(7), the density distribution, given by equation (jA.3|) . is 

n 7 = -|t 3 . (A.4) 

Far above their mass scales the massive particles are in thermal equilibrium. Therefore, 
the phase space distribution function is given by 

r = e"*, (A.5) 
where Ei = ^/\pi\ 2 + M 2 . So the equilibrium density distribution is 

< 9 =9i I Wf fr (A - 6) 

In terms of the dimensionless variables x = +m k anc j g _ M equation (jA.fij) can be 
rewritten as 

^3 r°° 



rpj POO 

n T = 9i — o / xe~ x Vx 2 — Z 2 dx 
^ Jz 



gi —Z 2 K 2 (Z) (A.7) 



'2tt 2 

where K 2 {Z) is modified Bessel function. As Z —>■ 0, z 2 K2(Z) — > 2. In this limit, the 
density distribution of massive particles resembles with mass less particles. Thus the 
approximation that at a temperature above their mass scale all particles are in thermally 
equilibrium is a valid assumption. As the Universe expands the temperature falls (Z 
increases). Therefore, the density distribution of all particles fall and is governed by the 
Boltzmann transport equations. 



A. 2 Boltzmann transport equations 

In this section, we give a pragmatic introduction to Boltzmann equations to study the 
evolution of the density of any particle species in the absence and presence of interactions. 

The expansion of the Universe dilutes the number densities of all types of particles 
even in the absence of interactions at a rate 

drii R 

~dT ~ 6 R ni 

= SHrii (A.8) 
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where R(t) is scale factor in Freedman Robertson and Walker (FRW) Universe and R is 
derivative with respect to time. H is Hubble expansion factor. Thus in the absence of any 
interaction the Boltzmann transport equation for the given particle species % of density 
rii is 

^ + 3Hn t = 0. (A.9) 

Now we scale out the effect of the expansion of the Universe by considering the evolution 
of the number of particles in a comoving volume. This can be done by dividing the number 
density of the particle species i with its entropy density, i.e. 

Y t = ^. (A.10) 

s 

Using the conservation of entropy per comoving volume (si? 3 = constant), equation (|A.10|) 
can be written as 

^ + 3Hm = sY = 0. (A.ll) 

As the Universe expands the momentum pi of the particle species i falls as 1/R and 
thus also the temperature T. Under rescaling the momenta of massless particles remain 
unchanged. So they keep themselves in equilibrium with the thermal plasma. Above the 
mass scale of any massive particle it will behave as a massless one. Below its mass scale 
the interaction rate decreases in comparison to the Hubble expansion rate and hence it 
falls out of equilibrium because it needs several collision times to keep it in equilibrium 
with the thermal photons. The departure of the density of any species i from its thermal 
equilibrium value can be predicted by solving the Boltzmann transport equations. 

For simplicity we consider the decay of any massive species i to a set of particles Y. 
As a result the equation (|A.11I modifies to 



^7(* - Y) - ^(Y 



rk- n 



Y 



(A.12) 



where 



l{i - Y) = J dU t dU Y (2^6\ Pl -p Y )ft q \A{i - Y)\\ (A.13) 

Inequation (lA~T3|) . dU = ^-^s- If we neglect CP- violation then \A(i -> Y)\ 2 = \A(Y -> 
i)\ 2 . Using (|A.6|) the above equation (|A.12|) simplifies to 

^ + 3 ^ = _r D (n i -< <? ), (A. 14) 
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r "=2k/p|k (2 ^ 4fe -^ )M|2 ' (A - 15) 

Note that in the above simplification we have assumed ny = n eq and it is true because the 
decay products Y are massless till the later epochs of our interest. Substituting Z = Mi/T 
and Yi = rii/s m. equation (|A.12|) we get 

*"i £g ly yeq\ 

dz zH{zy 1 i } 

= -D(Y i -Y*' 1 ). (A.16) 

Considering the 2 «-> 2 scatterings involving the species % equation (|A.16|) can be extended 
to 

dY 

-^ = -(D + S)(Y-Yn, (A. 17) 

where S = T S /ZH . This is the final Boltzmann equation for the evolution of any species 
i due to its decay and scatterings. 



where we have used 
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Appendix B 
Homotopy theory 



In the following we discuss some of the fundamental definitions which are useful for our 
purpose. 

Path : A path / in a manifold M is defined as a continuous function f(t) of a real 
parameter t, so that each value of t in the interval < t < 1 corresponds to a point 
f(t) in the manifold M. If a path f(t) connects the points P and Q we have f(0) = P, 
/(l) = Q. If /(0) = /(l) = P we have a c/osed path. On the other hand, if /(0) ^ /(l) 
then the path is open. 

Inverse of a path / is written as f^ 1 and is defined as 

f-\t) = f(l-t), (B.l) 

so that it corresponds to the same path traversed in the opposite direction. 
Product of two paths / and g is written as h = fg and is given by 

h{t) = f(2t) for 0<t< 1/2 

h(t) = f(2t - 1) for 1/2 < t < 1. (B.2) 

Two paths f(t) and g(t), both starting at P and ending at Q(P ^ Q), are said to 
be homotopic to each other if f(t) is continuously deformed to g(t) and is defined by 
/ ~ g. More specifically we can define a function L(s,t) such that L(0,t) = fit) and 
L(l,t)=g(t). 

Having defined path we shall construct a group by introducing a class of paths ho- 
motopic to /, denoted by [/]. They must, of course, have the same end-points. These 
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homotopy classes may be multiplied, by denning the multiplication law 

[/][<?] = [M (B.3) 

It is easy to see that, this multiplication law defines a group, called the fundamental group 
or first homotopy group of the manifold M, and denoted by 71*1 (M). This requires the four 
properties that satisfied by a group. They are: 

1. Closure: If [f] E m{M) and [g] E vr^M), then adhering to fffOfl [f][g] E tt^M). 

2. Associativity: Since (fg)h ~ f(gh) we have ([/][5<])[/i] = [f}([g}[h]). 

3. Identity element: This is the class of paths [I] that can be shrunk to a point. 

4. Inverse: Since [/-!][/] = [I], so [Z]" 1 = 

The manifold M is said to be a simply connected if all the closed paths in that manifold 
can be shrunk to a point. On the other hand, if all paths can not be shrunk to a point 
then, the Manifold is path connected. 
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